IHES/P/10/16 



Statistical Properties of Cosmological Billiards 

Thibault Damour 1 ' 2 'Q and Orchidea Maria Lecian 1 ' 2 > 3, 

1 Institut des Hautes Etudes Scientifiques, 35, route de Chartres, 91440 Bures-sur-Yvette, France. 

2 ICR AN et, Pescara, Italy 
3 APC, UMR 7164 du CNRS, Umversite Pans 7, 
10, rue A. Domon et L. Duquet, 75205 Paris Cedex 13, France. 

Belinski, Khalatnikov and Lifshitz (BKL) pioneered the study of the statistical properties of the 
never-ending oscillatory behavior (among successive Kasner epochs) of the geometry near a space- 
like singularity. We show how the use of a "cosmological billiard" description allows one to refine 
and deepen the understanding of these statistical properties. Contrary to previous treatments, we 
do not quotient the dynamics by its discrete symmetry group (of order 6), thereby uncovering new 
phenomena, such as correlations between the successive billiard corners in which the oscillations 
take place. Starting from the general integral invariants of Hamiltonian systems, we show how 
to construct invariant measures for various projections of the cosmological-billiard dynamics. In 
particular, we exhibit, for the first time, a (non-normalizable) invariant measure on the "Kasner 
circle" which parametrizes the exponents of successive Kasner epochs. Finally, we discuss the 
relation between: (i) the unquotiented dynamics of the Bianchi IX (a, b, c or mixmaster) model; (ii) 
its quotienting by the group of permutations of (a, b, c); and (iii) the billiard dynamics that arose in 
recent studies suggesting the hidden presence of Kac-Moody symmetries in cosmological billiards. 

PACS numbers: 04.50.-h Higher-dimensional gravity and other theories of gravity; 05. 45. -a Nonlinear dy- 
namics and chaos; 98.80.Jk Mathematical and rclativistic aspects of cosmology 



* damourSSihcs.fr 
1 lecian@ihcs.fr 



2 



I. INTRODUCTION 

A remarkable achievement of theoretical cosmology has been the construction, by Belinski, Khalatnikov and Lifshitz 
(BKL), of a general solution to the 4-dimcnsional vacuum Einstein equations in the vicinity of a spacelike ("cosmologi- 
cal") singularity pj, 0, 0- They found that this solution exhibits a never-ending oscillatory behavior, with strong 
chaotic properties. They could describe in detail the statistical properties of this never-ending oscillatory behavior 
by approximating the Einstein field equations (near the singularity) by a system of ODE's for three variables a, b, c 
('anisotropic scale factors'), namely 

J2 i„ „ 

a 4 , (1.1a) 
b\ (1.1b) 

(1.1c) 



d 2 In a 
dr 2 


= (b 2 


~c 2 ) 2 


o d 2 In 6 
dr 2 


= (c 2 


-a 2 ) 2 


n d 2 lnc 
dr 2 


= (a 2 


-b 2 ) 2 



where dr = —dt/(abc), and by approximately reducing the continuous dynamics of a, 6, c to a sequence of discrete 
maps. The crucial discrete map introduced by BKL relates the 'Kasner exponents' p a , pi,, p c describing the (ap- 
proximately linear) r-evolutions of the three scale factors a, b, c during two successive 'epochs' (i.e. two successive 
segments of the dynamics (|1.1[) during which the influence of the right-hand side is negligible). More precisely, BKL, 
following 0, parametrize the three Kasner exponents p a , pb, p c (constrained to satisfy 1 = p a +pb +p c = Pa ~\~Pb +Pc) 
by means of one real parameter u, and show that the interval 1 < u < oo is in one to one correspondence with the 
unordered set {p a ,pt,,p c }. In terms of this parametrization, BKL showed that the discrete map describing the passage 
from one epoch to the next is 

if u > 2, v! = u - 1 (1.2a) 

while, if 1<m<2, v! = . (1.2b) 

u — 1 

They also defined an 'era' as being a set of successive epochs during which u evolves according to the simple law 
(ll.2aj) . This led them to realize that the 'chaotic' part of the discrete epoch dynamics (jl.2[) is essentially contained in 
the "Gauss iteration map" 

X n+ i = {l/x n } = l/x n - [l/x n ]. (1.3) 

Here, x (with < x < 1) denotes the fractional part of u during an era. Let us recall that the u parameter of all the 
epochs belonging to the n-th era can be written in two essentially equivalent ways, which depend on the precise way 
in which one defines an era as a collection of epochs. If one uses the definition of an era such that the corresponding 
u parameters are always larger than one, the n th era consists of k n epochs parametrized by values of u of the form 
k n + k n — 1 + x n , 1 + x n . The next era then starts by a value of u equal to fc n +i + x n+ \ = l/x n , so that 
k n+ i = [l/x n ], and x n+ i = {l/x n }. We shall refer to this definition of an era as the "standard" one (because it 
was adopted in the treatise of Landau and Lifshitz Q), or as the BKL u> i one. An alternative definition, introduced 
by BKL in Eq. (5.4) of Q (that we shall call the BKL tl>0 one) leads to considering that the n-th era starts with 
u = k n + x„ — 1, and ends with u = x„ < 1 (with the next era starting with u = k n+ \ — 1 + x n+ i). We shall see 
below that the latter, (alternative) definition of an era is more natural in the billiard picture, and this is the one we 
shall actually use in our work. The iteration of the Gauss map (|1.3[) leads to a statistical behavior of the successive 
values of x, with an asymptotic stationary probability distribution over the interval [0,1] Q, 0, 0: 

1 dx 

w{x)dx = ^2TTx- (L4) 

For general reviews of the works dealing with the BKL singularity , see [1] 3- Later studies have refined the 
description of the statistical properties of the chaotic BKL oscillations, notably by introducing and studying more 
complete discrete iteration maps (involving several real variables), and notably by the two-dimensional discrete map 

1, 1, M 

<+i = {1/4} = V4 - [1/41, (l-5a) 
^+i = l/([l/4] + z»)- (L5b) 



1 Here {y} denotes the fractional part, {y} = y — [y] (where [y] denotes the integer part of y) of the (positive) real number y. We recall 
that the Gauss map is at the basis of the expansion of a positive real number into a continued fraction m + l/("2 + 1/(^3 + ■■■))■ 
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The iteration of this two-dimensional map (of the unit square into itself) asymptotically leads to a statistical 
behavior for e [0, 1] x [0, 1] with probability distribution (9j. [Toj. [TTj 

w{x + ,x-)dx + dx~ = d * d f C _. 2 - (1.6) 

In2 (1 + x + x ) z 

Separately from the work of BKI0, and, with a different aim and motivation, Misncr realized that generic Bianchi 
IX homogeneous cosmological models have a "very complex singularity" [l3|. He described this complex dynamics 
by a Hamiltonian formalism, in terms of a "system point" ft = (/3 + ,/3„) bouncing against a system of "potential 
walls". A reformulation of this dynamics [l4| . 15 1 led to the simpler picture of a point moving on a Lobachevsky 



plane and reflecting upon fixed billiard-type cushions. This led Chitre (using earlier mathematical results by Hopf 
and Hedlung) to remark that the dynamics of the system point is ergodic and mixing, with unique invariant Liouville 
measure (restricted to a fixed energy shell) , 

fi L = 6(H(q,p) - E)d 2 qd 2 p ex j^y- (1.7) 

The description of cosmological singularities in terms of billiards in (higher dimensional) Lobachevsky (or Lorentzian) 
spaces has recently received a new impetus from the discovery that the billiard chambers corresponding to many 
interesting physical theories can be identified with the "Weyl chambers" of certain (infinite-dimensional) Lorentzian 
Kac-Moody algebras [l?]], [H|, 0- This has raised the conjecture that, hidden below the BKL "chaos", there lies 
a remarkable "Kac-Moody symmetry" , akin to the duality symmetries of supcrgravity and string theories [2C| , [2lj , 

M- 

Coming back to the cosmological singularities in (3 + l)-dimensional General Relativity, the problem of relating the 
statistical properties of the discrete BKL map, such as Eq. (|1.4j) or Eq. (|1.6|) . to the invariance of the Liouville 
measure (|1.7[) in the continuous billiard dynamics, a la Misner-Chitre has been considered in some detail by Kirillov 
and Montani [23j . These authors have shown, by an explicit calculation, that the Liouville measure hl, Eq. (|1.7p . 
(which is a three-form) could be formally rewritten as the product of the invariant measure of the discrete BKL-type 
map (|1.6p . namely the two-form 

dx + A dx~ . 

/i2 = (i+x+x-r ( } 

by a one-dimensional measure dp, measuring the proper (hyperbolic) length along the billiard motion on the 
Lobachevsky plane. One of the aims of the present work is to better understand the link between the two different 
invariant measures (|1.7|) . (|1.8[) . and the origin of these measures within the symplcctic structure of the (Lorentzian) 
billiard dynamics. Another aim will be to go beyond the symmetry quotienting which has been used in most previous 
studies of the statistical properties of cosmological billiards. Indeed, there is a basic trialitjH symmetry between the 
three BKL dynamical variables a, 6, c, and the discrete maps (|1.3[) or (|1.5[) arise only if one effectively quotients the 
phase-space dynamics by their symmetry. [An example of this quotienting is the fact that the parameter u, taken in 
the interval 1 < u < oo, parameterizes the unordered set of Kasncr exponents {p a ,pb,p c }-] Here, we shall instead 
consider the richer (continuous and discrete) billiard dynamics in the full, unquoticnted phase-space. As we shall 
see, this full dynamics contains new statistical features that do not appear in the traditionally considered quotientcd 
dynamics. Finally, we shall also compare and contrast the (unquotiented or quotientcd) BKL dynamics of the (diag- 
onal type-IX) a, b, c system with the billiard dynamics that naturally arose in the recent studies that uncovered the 
hidden presence of Kac-Moody- related structures in cosmological billiards ([l^l, [3, 0' Hot [U, H3)- Indeed, the 
billiard dynamics which are most closely connected to such hidden symmetries take place in the Weyl chambers of 
some Kac-Moody algebras. In the usual case of 4-dimension vacuum Einstein gravity, this Weyl chamber is what we 
shall call a 'small billiard' , obtained by quotienting the full (a, 6, c) billiard table by the six-fold permutation group 
of the three letters a, b, c. As we shall discuss, the billiard dynamics in this quotientcd (configuration space) billiard 
table is not equivalent to the quotienting of the billiard dynamics in the full table, though the further quotienting of 
this small billiard dynamics by modding out the action (in phase space) of the a — b — c permutation is equivalent 
to the (phase-space) quotienting of the full 'big billiard' dynamics. Our present paper will focus on the usual case of 
4-dimensional vacuum Einstein gravity. In a sequel paper, we will extend our results to higher-dimensional gravity 



models, using the generalized "cosmological billiard" approach [21 1. 



2 It seems that western physicists, and notably J. A. Wheeler who was in the audience, first heard about the BKL results from a seminar 
given by Isaak Khalatnikov at the Institut Henri Poincare, Paris, in January 1968; see |12| |. 

3 Actually, what matters is a six-fold symmetry corresponding to the permutation group of the set {a, b, c}. 
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The paper is organized as follows. In Section II, we introduce cosmological billiards; this leads, in particular, to con- 
trasting the 'big billiard' studied by Bclinski, Khalatnikov, Lifshitz and Misner, with the 'small billiard', connected 
with the Wcyl chamber of a Kac-Moody algebra as shown in Fig[2j In Section III, we discuss two conformal represen- 
tations of cosmological billiards, i.e. the disk model and the upper-half-plane model; in particular, we outline their 
representations of epochs and eras. In Section IV, we define integral invariants for general dynamical systems, which 
allow us to find an invariant measure for the BKL discrete map (full details about integral invariants in Hamiltonian 
systems are given in the Appendix). In Section V, we analyze the big billiard, describe its dynamics as a "hopscotch 
game" at different levels, and define the corresponding maps. In Section VI, we exploit the symmetries of the big bil- 
liard to define a symmetry-quotiented map. In Section VII, we study the main properties of this symmetry-quotiented 
map. In Section VIII, periodic phenomena in cosmological billiards are considered, and some differences between the 
complete billiard and the symmetry-quotiented billiard are outlined. In Section IX, the small billiard is introduced: 
its features are investigated, and its equivalence with the big billiard is discussed. Brief concluding remarks end the 
paper. 



II. REMINDERS AND TECHNICAL PRELIMINARIES 



Let us start by defi ning our notation and recalling some basic facts about "cosmological billiards". [We mainly 
follow the notation of |2l|.] In order to describe the evolution of a general inhomogeneous space-time metric near a 
space-like singularity, it is convenient to use "pseudo-Gaussian " coordinates, with vanishing "shift" N l = 0, but with 
some convenient choice of the "lapse" N: 

ds 2 = - (N(x°,x l )dx ) 2 +g lJ (x ,x k )dx z dx : >. (2.1) 

Here x° denotes the coordinate time associated to any particular way of choosing the value of the lapse function 
N. The indices i, j = l,...,d denote the various spatial dimensions. In the present work we shall consider the case 
d = 3, but many of the general technical results recalled in this section are valid for any value of the space dimension 
d. There are two useful choices of the lapse N for exhibiting the "billiard nature" of the dynamics of gij near the 
singularity. The choice 

N = y/g, (2.2) 

where g denotes the determinant of the spatial metric g^, corresponds to using as coordinate time x° the parameter 
r introduced by BKL, i.e. the quantity 

dr = ~% (2-3) 

where dt = Ndx° denotes the (local) proper time. [A minus sign is introduced in (|2.3|) so that the cosmological 
singularity conventionally located at t — s- + ("big bang") is approached as t— s-+oo with respect to the r-coordinate 
time.] 

The choice of (|2.2j) . (|2.3j) leads to an asymptotic description of the gravitational dynamics in terms of a "Lorentzian 
billiard in /3-space" . More precisely, one first performs an Iwasawa decomposition of the spatial metric, i.e. one 
(locally) replaces the d(d + l)/2 functions gij(x°,x k ) by the d functions P a (x°,x k ), together with the d(d — l)/2 
functions Af a i (x°,x k ) parametrizing an upper triangular matrix with l's on the diagonal, according to 

d 

a=l 

In the near-singularity (or BKL) limit (t —> + or r — > +oo, or ^ a f3 a — > +oo), one finds that the upper triangular 
matrix has a limit [21| and that the only parts of the metric which have a "chaotic behavior" are the "diagonal 
degrees of freedom" parametrized by the d functions f3 a (x°, x k ). Then one finds that, at each point of space, the /3 a 's 
asymptotically follow a Lorentzian billiard dynamics: namely, the /3 a (r)'s undergo a succession of constant- velocity 
straight-line flights interrupted by collisions (and reflections) on some hyperplanes in /3-space. 
The free-flight dynamics of the (3 particle between wall collisions is described by the free action 



a, 6—1 



(2.5) 
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submitted to the constraint 

E^!=0- (2-6) 



a,b=l 

Here, the /3-space metric G a b is defined by 



da 2 = Gabd(3 a df3 b = £ (dft a f - d 0°) • ( 2J ) 

a,&=l a \ a / 

The metric G a & endows the d-dimensional space of the ft's with a Lorentzian structure (signature h +...+). 

[In the case of 3 + 1 dimensional General Relativity, the 3-dimensional ft space metric has signature — + +. 
Note, however, that the coordinates ft 1 , ft 2 , ft 3 are not of the canonical Lorentzian form. Indeed one has da 2 = 
—2 (dft 1 dft 2 + dft 2 dft 3 + dft 3 dft 1 ).] Note that the constraint (|2.6[) means that, between collisions, the ft particle goes 
"with the velocity of light" (in the sense of the Lorentzian structure G a b in ft space.) In other words, the free-flight 
dynamics deduced from the action (|2.5|) . namely 

^=0^ ft a (r)=ftZ + v a T, (2.8) 

is restricted by the quadratic constraint 

= £ G ab v a v» = J>") 2 - (j2 v a ) ■ (2-9) 



.6=1 



The free-flight dynamics (|2.8[) in ft space (and in the r parametrization) corresponds, in the BKL language, to a 
"Kasncr epoch" (between two successive wall collisions). The usually considered Kasner exponents p a , a = l,..,d, 
corresponding to a Gaussian (or synchronous) gauge (i.e. N = 1), are related to the d-dimensional velocity vector v a 
in ft space (and pseudo- Gaussian r gauge, N = y/g) via 

Pa = (2.10) 



Note that while v a satisfies the unique quadratic constraint (|2.9[) (proportional to the combination J2 a (Pa) 
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a p a ) 2 ), the Kasner parameters p a satisfy the two well-known constraints 

£p» = 1 = Ep2- (2-11) 

a a 

The free flight dynamics (|2.7p is only valid if the 'point' ft is sufficiently far from certain (Lorentzian) 'wall hyperplanes' 
in /3-space. The equations of these wall hyperplanes depend on the field content of the theory that one considers, (e.g. 
Einstein-Maxwell versus Einstein, etc.). They are of the general form 



G8) = £«^8» = 0. (2.12) 



More precisely, the dynamics of the ft particle is given by a Hamiltonian H of the form H = Hq + V, where Hq is a 
free kinetic term describing (in Hamiltonian form) the free-flight part of the dynamics, namely 

1 d 

Ho = ^J2 Gah ™, (2.13) 

a,b=l 

where G ab is the inverse of the covariant metric tensor G a b (in ft space), introduced in (|2.7[) . Its components in the 
ft a (Iwasawa-related) coordinates are explicitly given by 

J2 Gab ^»=J2M 2 -^i(l2^) ■ ( 2 - 14 ) 

a.b a \ a / 
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As for the potential V(/3) in the Hamiltonian H = Hq + V, it is a sum of "Toda-like" , i.e. exponential, terms: 
~ c a ex P( — 2w A (f3)). As recalled in the Appendix, in the near-singularity limit the potential V(f3) can be 
replaced by its sharp-wall limit Voo(fl) = J2a @oo(— 2w a ([3)) where Q(x) := if x < and @(x) := +oo if x > 0. 
For a generic inhomogeneous metric, the set of linear wall forms w A {[3) always includes curvature (or gravitational) 
walls {u!? abc \ (/3)) and symmetry (or centrifugal) walls (w^ b ). They are explicitly defined by 

wl abc) {P)^p a -p b -p c + Y,P e ( fe ^ c )> ( 2 - 15 ) 



«f o6) = (3 b - /3 Q (o<6). (2.16) 

Beware of the fact that the indices with parentheses appearing on the left-hand sides (l.h.s.) of these definitions 
should be considered as labels (like the label A in Eq. (|2.12|l ). and not as /3-space tensor indices. E.g. in the linear 
form of the /3's U!? abc -,(j3) = ^2 e w^ abc ^ e f3 e only the summed-over index e must be considered as a tensor index. In 

addition, note that the index e on w 9 r abc \ e is covariant, while the index e on /3 e is contravariant. This means that, 
when computing the Lorentzian scalar product between two wall forms, w A {[3) = ^2 a w A [3 a and w B ((3) = ^ a w B (3 a , 
one should use the contravariant /3-space metric G ab : 



w A -w B 



J2G ab w A w B . (2.17) 



a . b 



Among all possible walls entering the Hamiltonian, only the subset of "leading" walls (those not "hidden behind" 
another wall) should be retained to define the /3-space billiard defining the asymptotic BKL-like dynamics. Indeed, 
the billiard chamber is defined as the intersection of the positive sides of the set of wall hyperplanes, i.e. the 
domain where all the linear forms w A (f3) are positive. For instance, in the case of three spatial dimensions, there 
are 6 gravitational walls, and 3 symmetry ones. However, among these, some walls are "subleading" in that they 
are always behind some other walls. E.g. the symmetry wall w^ 3 ^ = /3 3 — /3 1 can be identically expressed as 

w (i3) = /? 3 — /? 2 + /3 2 — /3 1 = w? 12 \(/3) + w^ 23 ^(f3). Therefore the inequality w? 13 s(/3) > is a consequence of the 
two inequalities w? 12 s {(i) > and w^ 23 ^ (/?) > 0, meaning that the wall w^ 13 ^ (/3) is behind the two walls w^ 12 ^ ((3) and 
w (23) (P) anc ^ ^ s therefore subleading. Similarly, one finds that among the gravitational walls w 9 ^ abc ^ the ones where 
the first label is equal either to b or c (b ^ c), i.e. (J, C (P) = —f3 c + ^2 e P e are always subleading. In d = 3, this leaves 
only 3 a priori leading gravitational walls 

wf 123) (p)=2p 1 , ^ 231) (/3)=2/3 2 , wf 312) (p) =2/3 3 . (2.18) 

Moreover, the same reasoning which allows one to conclude that the symmetry wall w^ 13 ^(/3) is 'behind' the two other 
symmetry walls w^ 12 ^ (j3) and w^ 23 ^ (f3) shows that the gravitational walls w^ 231 ^ (f3) and w^ 312 ^ ((3) are both 'behind' the 
combination of the walls {w^ 12 y w^ 23 y w 9 (i 23 y!- Therefore, for a generic inhomogeneous metric the asymptotic billiard 
chamber in /3 space is defined by the following three independent inequalities 



W 



f 12) (^)=^ 2 -^>0, (2.19a) 
f 23) (/3)^ 3 -/3 2 >0, (2.19b) 
f 123) (/3)^2/3 1 >0. (2.19c) 



Note that the boundary of this billiard chamber is made of two (portions of) symmetry walls, and one (portion of) 
gravitational wall. The occurrence of the symmetry walls here comes from terms in the Hamiltonian associated with 
the kinetic energy of the off-diagonal components, , in the Iwasawa decomposition (|2.4[) of the metric. Note that 
an alternative way of seeing the 'constraining' effect of the off-diagonal components of the metric consists of using, 
instead of an Iwasawa decomposition, a Gauss decomposition of the spatial metric: gij = ^2 a e~ 2/3 R a i R a J 5 with 
R a i being a rotation matrix, parametrized by three Euler angles. Such a Gauss decomposition was introduced by 



Belinski, Khalatnikhov and Ryan in |2J] . As shown there it entails the presence of centrifugal walls which are simply 
related to the (Iwasawa) exponential symmetry walls exp(— 2w^ ab ^) via V^° & " tnf tx [sinh(w^ ab ^)}~ 2 . 
In the special case of a homogeneous vacuum model of Bianchi-typc IX, it is possible to restrict oneself (without loss 
of generality) to considering a metric gij9 l 6 3 which is diagonal in a co-frame 9 l — e l m dx m of left-invariant one-forms. 
In that case, the kinetic energy terms associated to the off-diagonal components of gij vanish, so that the symmetry 
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walls do not appear. As a consequence, the billiard chamber for the special diagonal Bianchi-IX case is defined by 
the three leading gravitational walls (|2.18[) . i.e. by the three inequalities 

(2.20a) 
(2.20b) 
(2.20c) 



g 

^(123) 


(/?) 


= 2/3 1 


>o, 
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W) 


= 2/3 2 


>o, 
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W (312) 


W) 


= 2/3 3 


> o. 



Note that these three billiard walls correspond to the leading terms that appear on the right-hand side (r.h.s.) of the 
BKL a, b, c system (| 1 . 1|) when using the exponential parametrization a — e~ a , b = , c ~ e~ 7 . Indeed, in terms of 
these variables, the a, b, c system fll.ljl reads 



d 2 a 
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e " 


-ia 




dr 2 


~~ 2 








d 2 (3 
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-4/3 


- e" 47 


dT 2 


- 2 








d 2 -f 
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e~ 
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2e 



-2(/3+7) 



the terms oc e 



-4a 



-4/3 



-47 



;-4/3 + 2e -2(/3+a) 



exactly correspond to the three wall forms (|2.20[) . i.e. e 4l3± , 



-4/3 2 



-4/3 3 



(2.21a) 
(2.21b) 
(2.21c) 
They appear 



even more clearly in the Hamiltonian constraint of the a, b, c system which has the form H = ^G a i,$ a $ + V((3) = 0, 
i.e., explicitly, 



H = 



da dj3 da dj d/3 dj 



dr dr dr dr dr dr 



,-4/3 



; -47 _ 2e -2« e -2/3 _ 2e -2« e -27 



2e- 2p e^] 



0. 



(2.22) 



It is easily seen that the three symmetry walls uij? 12 )(/3), w^ 3 ^(/3), w^ 31 ^(f3) partition the (i space chamber (|2.20|) 
into six sub-chambers which are all congruent (with respect to the Lorentzian geometry of /3-space) to the billiard 
chamber (|2.19[) corresponding to a generic inhomogeneous (and generically non-diagonal) metric. In view of this, and 
for brevity, we shall refer, in the following, to the diagonal Bianchi-IX chamber (|2.20[) as being the big billiard by 
contrast to the small billiard (|2.19p associated to a generic inhomogeneous metric. Note that, in the introduction, we 
referred to the big billiard either as 'the full a, 6, c billiard' or as the 'unquotiented a, b, c billiard'. In addition, note 
that the (six times smaller) billiard table of the small billiard is obtained by quotienting the big billiard table by the 
permutation group of {(3 1 , /3 2 , /3 3 }. 



A. Hyperbolic billiards 

So far we have recalled how the use of the time gauge (|2.2[) , (|2.3p leads to a description of the asymptotic dynamics 
of the metric, near a space-like singularity, in terms of a billiard motion in an auxiliary d-dimensional Lorentzian 
space parametrized by the 'logarithmic scale factors' f3 a . [Note that, in the diagonal Bianchi-IX case the BKL scale 
factors a, b, c are related to the /3's via a = exp(— /3 1 ), b = exp(— /3 2 ), c = exp(— j3 3 )]. A convenient reformulation of 
this Lorentzian billiard consists of decomposing the motion in /?-space into radial, p, and angular, r y a , parts. Here, the 
terms "radial" and "angular" refer to Lorentzian analogs of the usual Euclidean decomposition of a position vector x 
as x = rn, with r = (x 2 )^ and n being a unit vector. Namely, one decomposes the (time-like) 'position vector' f3 a in 
Lorentzian space as 

P a ^ P1 a ; P ^ (-G ab f3 a f3 b ) 1/2 . (2.23) 
When using such a decomposition, it is convenient to redefine the time gauge, and to replace the condition (|2.2p by 

N = p 2 J-g. (2.24) 

One then finds that the radial motion asymptotically decouples from the angular one and leads to a uniform motion 
of the logarithmic variable A = lnp w.r.t. the coordinate time, say T, associated to the new gauge (|2.24[) , i.e. 

,„ dr dt 

dT=-— = -^—. 2.25 
P 2 P 2 V9 
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Hi 



7=0, 



FIG. 1. The hyperbolic billiard on the unit hyperboloid T-L2- The 'big billiard' chamber defined by the three gravitational walls 
is sketched. 



As for the "angular motion" 7 a (T) on the (future) unit (d— l)-dimensional hyperboloid, say 

U d -i : G ab7 a 7 h = -1, (2.26) 

it is found to be asymptotically described by a hyperbolic billiard, i.e. by a succession of constant-velocity (in T 
time, Eq. (|2.25p ) geodesic flights on the unit hyperboloid Wd-i (|2.26p . interrupted by collisions (and reflections) on 
hyperbolic walls located on some geodetic hyperplane in 7-space. The 7-space chamber within which this angular 
billiard dynamics takes place is simply the projection (seen from the origin) of the corresponding /3-space billiard 
chamber onto the unit hyperboloid Hd-i, Eq. (|2.26[) . For instance, in the case of D = 4 vacuum Einstein gravity that 
we shall focus on in this paper, we end up with a non-compact, but finite volume billiard chamber on a 2-dimcnsional 
hyperboloid H2 bounded by three geodetic lines. In the case of a generic, inhomogencous metric, the billiard chamber 
is defined by the H2 projection of the inequalities (|2.19p . i.e. 

'small billiard' : 7 2 - 7 1 > 0; 7 3 - 7 2 > 0; 2 7 x > 0. (2.27) 

On the other hand, in the special diagonal Bianchi-IX case, the billiard chamber is the H2 projection of (|2.20p . i.e. 

'big billiard' : 27 1 > 0; 2 7 2 > 0; 2 7 3 > 0. (2.28) 

The link between the projected 'big billiard' chamber on and the corresponding big Lorcntzian billiard in f3 space 
is sketched in Fig. Q] 

We have defined here the billiards on in terms of the three components of the unit Lorentzian vector 7° (satisfying 
Gabl a l b = -2 ( 7 *7 2 +7 2 7 3 + 7 3 7 1 ) = -1, see Eq. (|2.7j) with d = 3). This is a hyperbolic analog of defining a 
billiard chamber on the unit sphere S2 by writing three linear inequalities iua(h) > in the three components n , 
n 2 , n 3 of a unit Euclidean vector (satisfying (n 1 ) 2 + (n 2 ) 2 + (n 3 ) 2 = 1). For many purposes, the knowledge of the 
linear forms defining the billiard walls (e.g. u>j? 12 ) = 7 2 — 7 1 for the first wall of the "small billiard" (|2 . 2T[) ) is 

all that is needed to compute most quantities of physical interest. For instance, the (hyperbolic-geometry) angle 9ab 
between the walls 1x^(7) = WAal a and 103(7) = WBal a is given by 

cos ab = — ■ (2.29) 

V"i • WAy/WB ■ WB 

where wa ■ Wb = G ab WA a WBb, Eq. (|2.1T[) . G ab denoting the contravariant metric, Eq. (|2.14|) . E.g. one easily checks 
that the three angles on the small billiard (|2.27p are 0, § and -|, while the three angles between the three sides of the 
big billiard (|2.28[) are 0, 0, 0. In addition, the law of reflection of the 7-space T-time velocity vector, say V a = dj a /dT, 
on a certain wall wa{j) is simply given (at the location of the collision, i.e. when wa{j) = 0) by 

7 - a = r _ 2 MM | (230) 

WA ■ WA 
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where wa(V) = WA a V a , and where w A = G ab WAb is the contravariant vector associated to the covariant components 
WAa entering the wall form wa{i) = WAal a - Eq. (|2.30[) relates /3-space vectors that are (at the location of the 
collision) all tangent to Hd-i- It is obtained by projecting the corresponding /3-space, T-time collision law [l7| 

=v «- 2 WA( ^, (2.31) 

W A ■ W A 

which, contrary to (|2.30p . involves time- independent vectors. 



However, for some purposes, it is convenient to use an explicit parametrization of the billiard dynamics on Hd-i = 
~H n by means of n = d — 1 intrinsic coordinates. This can be done in several ways. Let us first emphasize that the 
unit hyperboloid H n is a model of the n-dimensional Lobachevsky space (i.e. it is diffeomorphic to ffi" and has a 
constant sectional curvature —1). As shown long ago by Beltrami (see, e.g., the textbook [2|| and the review paper 
j2t|). Tin admits several useful representations in an n-dimensional Euclidean space. Among these, the conformal 
representations are useful because, as they preserve angles, they allow one to express the reflection law (|2.30p as 
a usual (locally Euclidean) reflection law of the (local) velocity vector on the wall wa- The two main conformal 
representations are: 

(i) the 'ball model', say B n , which represents T-L n as a unit ball x 2 < 1 in n-dimensional Euclidean space x G R n , with 
metric ds 2 = 4dx 2 /(l — x 2 ) 2 ; and 

(ii) the 'upper half-space model', or 'Poincare model', say V n , which represents H n as the half-space v > 0, u = 
(w 1 ,...,?/ 1 - 1 ) G R"" 1 , with metric 

ds 2 = d " 2 + dv2 . (2.32) 

The 'ball' conformal representation can be geometrically realized within the (n + l)-dimensional Lorentzian /3-space by 
stereographically projecting (from the 'South Pole' 7s, i.e. a center of projection located on the past unit hyperboloid) 
the future unit hyperboloid (|2.26p onto a (n-dimensional) hyperplane passing through the origin in /3-space. 
The Poincare model can be obtained by a suitable geometric inversion of the ball model. In both models, the geodesies 
of % n become Euclidean circles orthogonal to the boundary (the boundary being a unit (n— l)-dimcnsional sphere for 
B n and the plane v = for V n ), while walls, i.e. geodetic hyperplanes, become (n— l)-dimcnsional spheres orthogonal 
to the boundary. Note that the 'boundary' corresponds to the 'absolute' of H n , i.e. its domain at infinity (which 
corresponds to the future null cone, G a b/3 a /3 b = when replacing the 7 a 's with projective /3 a coordinates). In addition 
to these Euclidean space representations (which naturally define n coordinates, (x , x n ), or (u , u n ~ , v), on 
Hn), it might also be useful to coordinatize T~i n by means of the hyperbolic analog of the polar coordinates on a 
sphere. E.g. in the case n = 2, one can represent the metric on H2 as ds 2 — d0 2 + sinh 2 dd(j> 2 . 



III. CONFORMAL REPRESENTATION OF THE D=4 COSMOLOGICAL BILLIARDS 



In space-time dimensions D = d + 1 = 4, the walls reduce to (geodesic) lines on the 2-dimcnsional Lobachevsky 
plane %2- As recalled in the previous section, one can consider two different pure gravity billiards, in D = d + 1 = 4: 
(i) the big billiard (|2.28p (defined by an 'ideal triangle' on H2, i.e. a triangle whose three sides meet at infinity 
with pairwise vanishing angles) or (ii) the small billiard (|2.27p (which has angles 0, ^ and ^ and only one vertex at 
infinity). The most symmetric representation of the big billiard (which manifestly respects the symmetry group, of 
order 3! = 6, of the inequalities (|2.28|l ) is obtained by using a disk model centered at the point 7 1 = j 2 = 7 3 . See 
Fig. [2] which also exhibits the small billiard (|2.27l) . By the Gauss-Bonnet theorem, (A + B + C — it = J KdS) the 
(hyperbolic) area of the billiard is equal to 7r, while that of the small billiard is tt/6 (consistently with the fact that 
there are six congruent copies of the small billiard within the big one) . By using a Euclidean geometric inversion with 
respect to the 'cusp' (i.e. the vertex on the absolute) of the small billiard, one obtains a Poincare model of the billiard 
in which that cusp is represented by the point at infinity of the upper half plane (v = 00). In this representation (see 
Fig. [3]) the geodesies 7 1 = 0, 7 2 = and 7 2 — 7 1 = are all represented as vertical straight lines. 

As the explicit form of the transformations relating the original gravitational variables /3 a successively to 7°, and 
to its images in the ball and Poincare models, tend to be unwieldy and not very illuminating, let us sketch how the 
form of the final results can be obtained essentially without calculations, by using various geometric considerations. 
In any Poincare representation (say of coordinates it, v) a general wall = wa(j) = wai^ 1 + WA2I 2 + WA3J 3 must be 
(projectively) equivalent to the equation of a circle. Therefore each 7° must be of the form 7 a (u, v) = A(w, v)S a (u, v), 
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FIG. 2. The disk model B2 of the hyperbolic billiard. Both the big billiard (with walls a, b, c) and the small billiard (with 
walls G, B, R) are sketched. The 6 fundamental Kasner intervals are indicated on the boundary of the disk, which is identified 
with the Kasner circle. 




- 1 

FIG. 3. The Poincare model Vi of the hyperbolic billiard. Both the big billiard and the small billiards are sketched. 



with 

S a (u, v) = A a (u 2 + v 2 ) + B a u + C a v + D a , (3.1) 

for some constants A a , B a , C a , D a . Moreover, if we choose to put the cusp of the small billiard at infinity in the 
Poincare plane, <5 1 (u,w) = and 5 2 (u,v) = must be the equations of two vertical lines (see Fig. [3]). Therefore, S 1 
and S 2 must simply be of the form S 1 = B 1 u + D 1 , S 2 = B 2 u + D 2 . By contrast, S 3 (u, v) = must be the equation of 
a circle centered on the v = axis, i.e. 8 3 (u, v) = A 3 (u 2 + v 2 ) + B 3 u + D 3 , and cutting the v axis at the two points 
<5 1 (u) = and 8 2 (u) = (see Fig. [3}. In addition, the infinity of "H 2 projectively corresponds to the light cone in 
/3-space so that the infinity of the Poincare model (i.e. v = 0) must correspond to G a bS a S b = 0. All these conditions 
fix the expressions of S a (u,v) modulo an overall factor and modulo the parabolic subgroup of the symm etry group 
(SX(2,R)) of %2 leaving fixed the cusp: i.e. v! = au + b, v' = au. It was shown by Kirillov and Montani [231 that a 
particular choice of a and b leads to expressions for S a (u, v) which are nicely compatible with the u-parametrization 
of Kasner parameters which has been used by BKL g, [|. This choice leads to the expressions 

S^u^v)^-^ (3.2a) 
S 2 (u,v) = u + 1, (3.2b) 
S 3 (u,v) = u(u + 1) + u 2 , (3.2c) 
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which entail [in view of the quadratic constraint (|2.26j) ] 7 a (u,v) = S a (u, v)/(y/2 v), i.e. explicitly 



7 1 (u,u) = 
j 2 (u,v) = 



V2v' 
u + 1 

72V 

u(u + 1) + v 2 
V2v ' 



As shown in Fig. [31 in this normalization the gravitational wall 7 1 
small and the big billiard is located at u = 0; the symmetry wall 7 1 = 
u = —1/2; and the other vertical gravitational wall of the big billiard, 7 2 = (b wall), is at u = —1. On the other 



(3.3a) 
(3.3b) 
(3.3c) 

= (i.e. the o wall) which is common to the 
7 2 of the small billiard (B wall) is located at 



hand, the remaining walls (either 7 = or 7 



7 



0) are circles orthogonal to the v axis. 



An essential role will be played in the following by the images in the Poincare model of the 'Kasner epochs' of the 
cosmological billiard, i.e. the free nights between two successive wall collisions. In /3-space, these free-flight segments 
are described by uniform motion (in r-time) and in straight line, see Eq. (|2.7p . The (/3-space) 'velocity' of these free 
flights is described by the Lorentzian vector v a , submitted to the constraint of being null, Eq. (|2.9[) . The 'Kasner 
parameters' p a of each Kasner epoch are (projcctively) related to the components of the /3-space velocity v a by the 
relation (j2~TU)) . 

The disk-model (B2) projection of consecutive /3-space free flights is made of geodesic segments in £>2, i.e. arcs of 




FIG. 4. Kasner epochs of the big billiard in the disk model B2 of the hyperbolic billiard. 



circles orthogonal to the boundary circle of 82- Sec Fig. [U which represents the 'reflection' of these geodesies on the 
three gravitational walls of the big billiard. 

When represented in the Poincare model, instead of the disk model, each /3-space straight-line segment /3(t) = 
/3 + v a r (with n < t < r 2 ) gets mapped into a geodesic segment of the Poincare plane V2 , i.e. a segment of a circle 
orthogonal to the horizontal axis v = 0, say 

u = 2^ u+ + ~ 2^ U+ ~ cos ^' (3.4a) 
v = i I u + - u~ I sin0, (3.4b) 

with 0<6i<9<92<ir Here 9\ (which corresponds to 7*1 in the /3-space "upstairs" ) corresponds to the last collision, 
and 6> 2 (o r 2 ) to the next one. See Fig[SJ which represents the big billiard in the Poincare model V 2 . Here, we are 
considering oriented circles whose formal extension to the full interval < 9 < tt would start, when 9 = 0, at the 
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FIG. 5. Kasner epochs of the big billiard in the Poincare model Vi of the hyperbolic billiard. 



location u = u on the v axis, and end, when 9 = 7r, at the location u = u + on the v axis. [Note that the radius of 
the circle is \ \ u + — u~ |, while its center is located at u c = \{u + +u~), v c = 0.] 

For a given billiard table, the oriented pair of real parameters (u + ,u~) (exemplified in Fig. [5]) uniquely determines 
the (oriented) geodesic segment corresponding to some Kasner epoch. More precisely, it is easy to see geometrically 
that the 'end' parameter u + uniquely parameterizes the family of /3-space straight-line segments j3 a = (3q + v a r 
(without considering their r parametrization) that share a common (formal) asymptotic direction oc v a . In other 
words, u + uniquely parameterizes the three Kasner exponents of the considered Kasner epoch. The precise technical 
link between u + £ R and the three p a 's is obtained by relating successively p a : (i) to v a (see Eq. p.lOp 1 ); (ii) to /3 q (t) 
in the formal r — > +oo limit (via Eq. (|2.7|) ): (hi) to 7°(r) = /3 a /P in the same limit, and thereby (iv) to j a (u,v) in 
the corresponding (formal) limit u — > u + , v — > corresponding to the future endpoint of the circle in the Poincare 
model representing the formal extension of the Kasner epoch. In other words, 

Pa = — r = hm — ; , ^— = hm — r— — = hm — t- -. (o.o) 

Using the Poincare model expressions (|3.3|) , this leads to the result 

Pa =p B a KL (u + ) 1 (3.6) 

where the three functions p BKL (u) are the well-known w-parametrization of Kasner exponents introduced by Bclinski, 
Khalatnikhov and Lifshitz, namely 



u 2 + u + 1 ' 
u+ 1 
u 2 + u + 1 



In other words, as was found in Rcf . (23| , the BKL u-parameter can be interpreted as the location on the real axis 
of the future end point of the circle representing the considered Kasner flight in a suitably defined Poincare model of 

On the other hand, if we now consider the formal extension of the considered Kasner flight towards the 'past' (in 
r-time), there occurs a subtlety: the past end point u~ of the corresponding Poincare circle does not correspond 
(as one might naively expect) to the formal r — > — oo limit of /3 a (r) = /3q + v a r, but only to the finite past limit 
t -> t such that /3Kasncr( T ) = /?o + "t^ 1 " intersects the /3-space (future) light cone: G ab /3^ asncr (ro) / 3K asncr (T ) = 0, i.e. 
p{ T o) = 0. As this limit again corresponds to a point at infinity for 7Kasncr( r ) = /^Kasner i T )/p( T ) S 'rii, one finds, by 
changing the limits in Eq. (|3.5p that the past end point u = u~ , v = on the (oriented) Poincare circle parametrizes 
/?Kasner( T o) m a projective manner: 

asncr 

(to) 



as i.i ci: 

(to) 



Pa L (u-)- (3.8) 
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fcl 


< u < 1 


u = 1/w 


pi (it') = pi(u) 


P2(tt') = pz{u) 


P3(u') = P2(w) 




-1/2 < u < 


u' = — (1 + «)/« 


Pl(u') = P3(«) 


p 2 (u) = pi(«) 


p 3 (u') = p a (it) 


fe 3 


-1< u < -1/2 


u — —u/(u + 1) 


pi(tt') = P3(u) 


P2(u) — P2(u) 


P3(u') = Pi (it) 


A'4 


-2 < u < -1 


it' = -l/o + i) 


Pl(u') = P2(ti) 


Pz(m') = P3(ti) 


p 3 (u') =pi(w) 


A-.-, 


— oo < u < — 2 


It' = — Ji — 1 


Pl(ft') = Pa(u) 


p 2 (u') = Pi (u) 


P3(«') = P3(lt) 



TABLE I. Kasner transformation, fei, k$ mapping the indicated intervals of the u line ("Kasner circle") onto the fundamental 
interval 1 < u < +oo. 



As we shall need them below, let us note at this stage some features of the BKL u-paramctrization of Kasner exponents 
(|3.7j) . First, let us emphasize that the manifold of Kasner parameters, p a , restricted by the two constraints (|2.11|) . 
is topologically a (d — 2)-dimcnsional sphere [indeed Eqs. (|2.11|1 ~) represent the intersection of a (d — l)-dimensional 
sphere by an hyperplane]. In the case considered here where d = 3, this means that the three Kasner parameters 
pi, P2, P3 run over a topological circle. In fact, this 'Kasner circle' can be identified, in the disk representation of 
V.2, with the boundary of the unit disk, i.e. with the absolute of H 2 - See Fig. [5J In particular, the it-parameter in 
Eqs. (|3.7p should be considered as running on the extended real line R = R U {oo}. The extended line u G R is then 
naturally divided into the six permutations of the three letters (pi, p 2 , P3) which constitute the symmetry group of 
the two Kasner constraints (|2.11l) . As the latter symmetry group is generated by reflections in the symmetry walls, it 
is natural to divide the u line, i.e. the boundary of H2, by means of these symmetry walls: j 2 — 7 1 = 0, 7 3 — 7 2 = 0, 
7 3 — 7 1 = 0, that is, by considering the solutions of the equations pf KL (u) — p2 KL (u), or p^ KL {u) = p^ KL (u). or 
pf KL (u) = p^ KL (u). This leads to dividing the u-line into the intervals: 

(-oo,-2); (-2,-1); (-1,-5); (-^,0); (0,1); (l,+oo). (3.9) 

If, following BKL, we consider the interval 1 < u < +00 as a 'fundamental interval', over which the Kasner exponents 
are ordered as pf KL (u) < pf Ki (u) < p§ KL (u), the other five possible orderings of p\, p 2 , P3 will be obtained 
by applying to the variable u a transformation implementing a composition of geometric reflections accross some 
symmetry walls. It is well known that a geometric reflection acts on the complex variable z = u + iv of the Poincarc 
plane according to 

, az + b , _ , 

z' = — — , (3.10) 

cz + a 

with a, b, c, d G R and ad — be = +1. When acting on the boundary of the Poincare model (v = 0), and composing 
several reflections, this leads to transformations of the form 

u , = ± au + b a d-bc=l. (3.11) 
cu + d 

The explicit expression of the five transformations u' — fi(u), i = 1, ...,5 of the form (|3.10j) that map the first five 
intervals (|3.9|) into the last ('fundamental') one are given in Table HI In the following, we shall refer to these as 
"Kasner transformations". Note that each boundary Kasner transformation v! = fi(u) uniquely determines the way 
the corresponding combination of reflections acts on the interior of the Poincare model: if the sign in Eq. p. lip is 
+ (even isometry) it is z' = (az + b)/{cz + d), while if the sign in Eq. (|3.11|) is — (odd isomctry), it is given by Eq. 
p.lOp . If one adds the identity transformation, say ko (ko(u) = u), the six transformations {fco, k\, k§} constitute 
(under their composition) a realization of the permutation group of three objects (say, the three walls a, b, c). 



IV. REDUCED SYMPLECTIC FORM FOR BILLIARDS 



Let us briefly recall various ways in which one can define integral invariants for Hamiltonian systems. (For more 
details, see, e.g., [27[, [13 and the Appendix below) For a general (possibly time dependent) Hamiltonian dynamics, 
with (Hamiltonian) action (i = 1, ...,n) 

S = J dtlpitf - ff(q,p,t)] , (4.1) 
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one can use the Poincare-Cartan two-form 

Jpl = d 1 l A d P^ -dtA dH(q, p, t) (4.2) 

i 

to define a whole hierarchy of integral invariants of the (unparametrized or parametrized) Hamiltonian flow. Leaving a 
more general discussion to the Appendix, we shall consider here the case of a time-independent Hamiltonian H{q,p) 1 

and focus on invariants of the flow on a given (2n — l)-dimensional energy hypersurface, say £^ n , satisfying 
H(q,p) = E. For this case, the invariance of the energy-shell reduced Liouville measure ^ oc S(H(p 1 q) — E)dq 1 A 

dq 2 A ... A dq n A dp\ A dp 2 A ... A dp n is well-known. Less well-known is the construction of integral invariants based 
on the existence of a reduced symplectic 2-form 



(2) 



dq l A dpi 



(4.3) 



defined on the quotient Q^ n ^ = £^ n /J^h of the energy hypersurface H{q,p) = E by the unparametrized 
Hamiltonian flow. In other words, Q^ n 2 ^ = £ E 2n ^ /J~h is the (2n — 2)-dimensional space of unparametrized 

Hamiltonian motions on £g . [The reduced symplectic form (|4.3p is linked to the general theory of reduction 
of phase spaces with symmetry; it was used in Ref. [29j as a way to define a measure in cosmology] A concrete 
representation of the quotient space Q^ n 2 ' can be obtained by considering any "transverse section" Tu of £^ n , 

(2) 

i.e any "initial conditions" for the (unparametrized) Hamiltonian flow. The fact that "descends" to the quotient 
means that its restriction to any transverse section of Th is independent of the choice of section. 
As a concrete example of the reduced (symplectic) form (|4.3[) on the space of motions (or initial conditions), one 
can have in mind the symplectic form on the manifold of (unparametrized) straight lines in a Euclidean plane. A 
straight line L, i.e. x(s) (where s measures the length along the line), can be parametrized by two vectors, b, 
n, submitted to the two constraints: n 2 = 1, b • n = 0, namely x(s) = b + ns. In any Cartesian coordinate 
system, one can explicitly parametrize (b,n) by two real numbers: b = (— & sin a, b cos a), n = (cosa,sina), where 
b is the impact parameter between the origin and L, and a is the angle between the x axis and L. The reduced 
symplectic form on the 2-dimcnsional manifold of (unparametrized) straight lines L can be obtained by starting 
either from the unparametrized action Si = J \J dx 2 + dy 2 , or from the parametrized one S2 = J dt^(x 2 + y 2 ) 
submitted to a fixed energy constraint. For instance, S 2 corresponds to a 4-dimensional phase space {x, y,p x ,p v ) with 
Hamiltonian H = | (p 2 +Py)- Let us consider the 2-form (|4.3p reduced both by the energy constraint H =1/2 (i.e. 
Px+Py = 1) or Px = cosa, p y = sin a), and by restricting it, e.g., to the section = y(so) = bcosa + Sosina, i.e. 
so = — &cota. These two constraints reduce the phase space to a 2-dimcnsional one parametrized either by Xq,Px, 
where xq = [x] y —o = — fosina + sq cosa = —bj sina and p x = cosa, or by (b, a). The symplectic form reduced to this 
section yields 

w' 2 ] = dx A dp x = d ( — ) A rf(cos a) = db A da. (4.4) 

\ sma / 

We see that the impact parameter b and the angle a constitute canonical coordinates on the 2-dimensional manifold 
of straight lines. The reduced symplectic form (|4.4j) can also be written in a form that is manifestly invariant under 
the group of Euclidean symmetries, namely 

w red = rfb A n = db x A dn x + db y A dn y (4-5) 

with the algebraic constraints n 2 = 1, b • n = 0. It is also easily checked that any other section (e.g. x = instead 

of y = 0) yields the same reduced form. Note that, in the present example, the reduced phase space of Euclidean 

(2) 

straight lines is 2-dimcnsional so that the reduced symplectic form ur ^ furnishes directly a measure on the space of 
straight lines. 



We spent some time on this simple example because we shall be interested, in this paper, with the hyperbolic-plane 
generalization of this structure. More precisely, we shall find useful to consider the reduced symplectic form (|4.3p in the 
case where the Hamiltonian dynamics is that of a (radially-projected) billiard motion on 7-space, i.e. on the hyperbolic 
space T-l n . More precisely, we consider the case where n = 2, and use the Poincarc half-plane model. The corresponding 
phase-space is 4-dimcnsional, say (u,v,p u ,Pv) in the Poincare model, and the corresponding Hamiltonian reads 



H 7 (u,v,p u ,p v ) = -v 2 (pl+pl) + Vaodfav)). 



(4.6) 



15 



is simply 2-dimensional. Therefore, in that case, the general fact that the ambient symplectic form — da^ c can 

(2) 

"descend" onto the reduced symplectic form ujj red on the abstract quotient space Q 7 (and, therefore, on any section 



The energy surface £^ = {H 7 (u,v,p u ,p v ) = E 7 } is 3-dimcnsional. Finally, the quotient space Q 7 = Se^/J^h^ 

thi 
forr 

( 2) 

of the Hamiltonian flow) means that 0J^ Ted directly provides an integral invariant of arbitrary "snapshots" of the 

(2) 

billiard motion. In particular, if we take the snapshots corresponding to collisions on the wall, we conclude that uj^ led 
is invariant under the discrete billiard map T ■ Moreover, as we are considering here a case where the quotient space 

(2) 

Q 1 is two-dimensional, w 7r g d directly defines an invariant measure of the discrete billiard map T ■ [See the Appendix 

for a general discussion of the link between the measure on Q^ n ~ 2 ^ associated to u>j£) ed and the energy-shell Liouville 
measure.] 

(2) 

There are many ways to compute the reduced symplectic form u>^ led - Let us first note that it is the generalization 
of the measure discussed above, on the 2-dimensional manifold of straight lines in a Euclidean plane. Here, indeed, 

(2) 

k-'-yred * s a measure on the 2-dimensional manifold of geodesic lines in a Lobachcvsky plane. As in the Euclidean 
calculation above, Eq. (|4.4[) , 0J^ ted can be obtained by reducing the ambient symplectic form du A dp u + dv A dp v by 
two conditions: the energy-shell condition v 2 (pf l + p%) = 2E 7 , and a cross-section condition locally transverse to the 
Hamiltonian flow. It is easily seen that, if we consider for simplicity the energy shell E = ^ for a geodesic (i. e. that 
the geodesic motion on H2 proceeds with unit speed), the general geodesic line in £ E ^_i (i.e. a circle orthogonal to 
the measure boundary) can be parametrized as 

u = U - VcosO, v=\V\sinO, (4.7) 

sin 9 1 cos 9 cot 9 

p u = sign (V) = — , p v = = — -, (4.8) 

v V v \V\ 

where 9 grows, according to d9/dT = sin#, from to ir as T formally varies from —00 to +00. The two constants of 
integration (U, V) parametrize the two-dimensional manifold of geodesic lines on the Poincare half plane. They are 
related to the parameters u + , u~ used in Eq. (|3.4p above via 

U=^(u + +u~), V = ±(u+-u-). (4.9) 

(2) 

If, for instance, we use as cross-section to restrict ui\ = du A dp u + dv A dp v any v = const slice (i.e. = dv = 
sin9d\V\ + \V\ cos9d9) we get 



u^l = du A dp u =d{U-Vcos9)Ad[ -I = - "" y2 "' ■ (4.10) 



1 \ dUAdV 



Rewritten in terms of (u + ,u ), Eq. (|4.9p . this reads 

(2) _ du + A du 
(w+ — u~) 



<d = 27-T— -2- ( 4 - U ) 



The result (|4.10[) or f|4. 1 1[) is similar to the Euclidean measure (|4.4I) or (I4.5[) on the manifold of Euclidean straight lines. 
Analogously to the fact that the measure (|4.5[) was invariant under the group of Euclidean symmetries (translations, 
rotations and reflections), the measure (|4.10[) , (|4.11[) is invariant under the group of symmetries of the hyperbolic 
plane. This group is (R) x Z2, and it acts on the boundary of the Poincare model (i.e. on the parameters 
u + and ur) by transformations of the form (|3 . 1 1 [) . This group of transformations is generated by u'± = au± + 6, 
u' ± = —l/u± and u' ± = —u±. It is then easily seen that (|4.11j) is indeed invariant under each one of these generating 
transformations. Note that this also gives a direct proof that the reduced symplectic form (|4.11[) is invariant under 
the Hamiltonian flow of the billiard. Indeed, this flow is made of two types of evolutions: (i) a free-flight evolution 
during which u + and u~ do not vary, and (ii) collisions on the walls, during which the geodesic undergoes a hyperbolic 
reflection, i.e. a transformation of the type (|3.10p leading to Eq. (|3.11[) with a minus sign. 



V. HOPSCOTCH DYNAMICS OF THE BIG BILLIARD 



We shall start our investigation of the various possible cosmological billiards for pure gravity in d — 3 spatial 
dimensions by considering the 'big billiard table' delimited by the three gravitational walls appearing in diagonal 
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homogeneous Bianchi IX cosmological models, i.e. the walls w^ 123 ^(/3), w^ 231 ^(/3) and w^ 312 ^(/3), Eq. (|2.18[) . which 

respectively correspond to the terms a 4 = e _4a , b 4 = e _4/3 and c 4 = e~ 47 in the usual BKL representation [see Eqs. 
(|2.2ip . ()2.22|) ]. Using Eqs. (|2.18[) . which express 7° cx f3 a in terms of the coordinates (u,v) of the Poincare model, 
we see that, in this model, the 'a wall' (or /3 1 wall) is located along the vertical line u = 0; the '6 wall' ((3 2 ) along the 
vertical line u = — 1; while the 'c wall' (/3 3 ) is located along the circle u{u + 1) + v 2 = 0. The reflection laws of a point 
z = u + iv E V2 through these geodesies must have the form p,10[) . One can determine the values of the coefficients 
a, 6, c, d entering the transformation (|3.10[> by requiring that the transformation leaves point-wise fixed the circle 
through which one is 'reflecting' (in a hyperbolic- geometry sense) the point z = u + iv. Indeed, the condition z' = z, 
i.e. z = ~(az + b)/(cz + d) yields as a locus of fixed points the circle 

= czz + dz + az + b = c(u 2 + v 2 ) + (d + a)u + (d - a)v + b (5.1) 

which degenerates to a straight line when c = 0. 

By successively identifying the point-wise fixed circle (|5.1[) to the three diagonal Bianchi IX walls, one finds the 

following reflection laws (acting on the end points u , according to Eq. p. lip with a minus sign). For the 'a wall' 
(u = 0), -> A^) with 

A{u ± ) = -u ± - (5.2) 

for the l b wall' (u = -1), ->• B^) with 

B^) = -u* - 2; (5.3) 
and for the 'c wall' (u(u + 1) + v 2 = 0), — > C^m*) with 

Note that all these reflection laws act diagonally (i.e. separately) on u + and u~ . A billiard motion in the presently 
considered big billiard is a succession of geodesic flights (or 'Kasner epochs') connecting two different walls. For 
instance (as illustrated in Fig. [5]) 

...c ^fe^a— s-a— >c— (5.5) 

In this work, we shall define a 'Kasner era' as a set of Kasner epochs joining the same two walls, with the condition 
that the epochs preceding and following the considered Kasner era involve the third wall. The length of an era is 
defined as the number of epochs (i.e. geodesic flights) it contains. For instance, in the sequence (|5.5[) there is an era 
of length 3 between the walls a and b, namely 

E 3 (b,a) : b -> a -> b -> a. (5.6) 

Here, we have introduced the notation Ek(x, y) for an era of length k = 1, 2, 3, ... whose first free flight is from the 
wall x to the wall y (were x, y € {a, &, c}, x 7^ y). Note that (x, y) is an ordered pair as one should distinguish an era 
which starts on x and then goes to y, from an era which starts on y and then goes to x. Note also that the beginning 
of an era is defined by checking (as sketched in (|5.5[1 ) that the previous connecting flight started from a wall z 7^ x 
and y (similarly for the end of an era). Note that the above definition of a Kasner era corresponds to one of the two 
different^ definitions of an era considered by BKL. More precisely, it is the definition they consider in Eq. (5.4) of the 
review In this definition (called BKL tt> o in the Introduction of this paper), the BKL u parameter varies (using 
the notation of the Introduction) from k — 1 + x to x (with < x < 1 and k G K). [This contrasts with the other 
(more standard) definition of an era used by BKL, the BKL u> i in which the u parameter varies from k + x to 1 +x 
(so that u stays in the interval [1, +00])]. As noted in the second footnote on p. 753 of [3], the former (less standard) 
definition of an era (that we shall use here) is more natural when considering the dynamics of the variables In a, In 6, 
Inc. In terms of the billiard picture, this more natural character does correspond to the definition we gave above of 
collecting all the epochs joining the same two walls. By contrast, in the other BKL definition (k + x to 1 + x), the 
era corresponding, for instance, to the sequence (|5.5[) would consist of the three epochs c — > b — > a — > b, and what is 
in our definition the last epoch of the era (of the type b — > a) would be considered as the first epoch of the next era. 
Note that in both definitions the Kasner era has the same length: in our example a length 3; the last epoch c — >• b of 
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the preceding era having been added as a first epoch, in replacement of the last b — >• a in Eq. (|5.6[) . 

The shortest possible length of an era is k = 1, i.e. an era corresponding to only one epoch. E.g., as we shall see 
later, the simplest periodic big-billiard orbit (involving the golden ratio) proceeds along the equilateral (hyperbolic) 
triangle geodesically connecting the 'middles' of the three a, b, c walls, and is made of only one-epoch eras, say 

...c — >c— >a— > c..., (5.7) 

or the reverse. 



A. Hopscotch dynamics 

As just recalled, the dynamics of the big billiard is described as a sequence of eras Ek(x, y), where each era Ek{x, y) 
is made of k epochs, i.e. k 'arrows', x — > y — > x — > connecting the walls x and y (x ^ y,x,y € {a, 6, c}). Note 
that the last wall involved in the era E^{x,y) will be y if k is odd, and x if k is even. To describe mathematically 
the discrete big-billiard dynamics induced by the effect of the successive collisions, i.e. of the corresponding discrete 
'billiard maps' 7i a transforming the phase-space variables on a a; wall (just after the cc-collision) to the phase-space 
variables on the following y wall, it is useful to use the Poincare plane variables. [Note, however, that though 
the Poincare-plane variables are algebraically more convenient, it is generally more enlightening to geometrically 
visualize the billiard dynamics on the disk model. Sec Fig. SJ] As recalled above, in the Poincare plane each epoch 
trajectory (i.e. each geodesic segment connecting two successive sides) is uniquely parametrized by an ordered pair 
(w~,u + ), where u~ (respectively u + ) is the end point (resp. starting point) on the v = axis (or 'absolute') of 
the corresponding, extended geodesic. In terms of these variables, the billiard dynamics induces a discrete map T 
transforming a point in the (u~ ,u + ) plane (describing some epoch) into another point (u ~ ,u + ) (describing the next 
epoch): (it ~,u + ) = T(u~ ,u + ). Looking at Fig. [5j it is clear that the knowledge of (it~,M + ), i.e. the knowledge of 
the initial epoch, uniquely determines the wall on which it will next collide, and therefore uniquely determines the 
explicit form of the transformation T, among the three possible explicit forms A, B, C, listed in Eqs. (|5.2I) . (|5.3I) and 
(|5.4[) . It is also clear that the transformation T is one-to-one because its inverse T~ 1 is defined by "reversing the 
time evolution", i.e. exchanging the roles of u + and u~ . As for the iteration of T, T°ToT° -, it corresponds to 
composing a sequence of (u~,u + ) transformations [among Eqs. (|5.2|) . (|5.3|) and (|5.4[) ] corresponding to a sequence of 
wall collisions. E.g. the sequence (|5.5|) will correspond to successively composing the actions of 

...C^B^A^B^ A^C^ ... (5.8) 

on the (u _ ,u + ) plane, i.e. the combined map (in reverse order) 

...Co AoBo AoBoC... (5.9) 

Note also that the composition of maps (|5.9[) can also be expressed as the corresponding matrix product of the matrices 




(with ad — be = 1) corresponding [via u' ± = —(au± + b)/{cu± + d) = (—au± — b)/(cu± + d)] to the fractional linear 
transformations (|5 . 2[) . i.e. 

...C.A.B.A.B.C... (5.10) 

where 




(5.11) 



and where the dots in Eq. (|5.10[) denote the ordinary matrix product. 

Summarizing so far: the representation of the big-billiard dynamics in the (u~,u + ) plane is the following: (i) 
during each epoch (i.e. free flight) the reduced phase-space point (u~ , u + ) stays fixed; (ii) the effect of each collision 
on the wall a, b or c consists in transforming the phase-space point (it~, u + ) into a new point (u'~ , u' + ) = l~(u~ , u + ), 
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where the explicit expression of (it' , u' + ) is uniquely defined by the initial phase-space point (u , u + j3, and is either 
of the form (A(u~), A(u+)), (B(u~), B{u + )) or (C(u~), C{u+)) [with A^), B^), C^) given by Eqs. (IQl . 
(|5.3p and (|5.4j) . respectively], where the choice between A, B or C is determined by the a, b, c wall that is next 
crossed by the oriented geodesic defined by (u~,u + ). In other words, we can think of the (u~,u + ) plane as a big 
'hopscotch courto on which the representative phase-space point (u~ , u + ) jumps around, in a deterministic manner, 
(u~ , u + ) — > T{u~ , u + ) -> 7° T(u~ , u + ) — > ... according to a sequence of 'jumps' whose concrete form is of the type 
(|5.8[) . These jumps act diagonally, i.e. in the same way on u + and u~. 

As we shall discuss below, though each "collision transformation" A, B or C acts on u + (respectively, on u~) 
independently of u~ (respectively, of tt + ), one needs to keep track of the successive values of the pairs (u~,u + ) to 
determine the entire (two-sided) sequence of maps, such as Eq. (|5.9[) . corresponding to the billiard dynamics. The 
big-billiard hopscotch dynamics T just defined differs from the usually discussed BKL dynamics in several respects. 
Indeed, in order to simplify their discussion, and go to the essence of the Bianchi IX dynamics, Bclinski, Khalatnikhov 
and Lifshitz did not keep track of the order of the Kasner exponents during an era (i.e., in their notation, whether 
{PhPm) is (pi,P2) or (p2,Pi) in an era of oscillations between the a and b walls). Moreover, BKL further simplified 
their discussion by using the 6-fold permutation symmetry among (a,b,c), so that they also did not keep track of 
which unordered pair {a, b}, {b, c} or {c, a} an era referred to, nor of the ordering of the first pair in a given era. By 
contrast, our description above explicitly keeps track of both ordering and labeling issues. Below, wc shall discuss 
how one can 'quotient', in a precise manner, the more complete big-billiard dynamics down to the usual BKL discrete 
dynamics. 

As an example of fuller description in the big-billiard representation, note that the particular lcngth-3 era E%(b, a), 
Eq. (|5.6[) . corresponds to applying, successively, to the values {u^ ba ,Up ba ) parametrizing the first epoch of £3(6, a) 
(i.e. the first arrow in (|5.6p ) the transformations A,B, and then A. This yields successively 

«± = ^(«!j = -4i«> ( 5 - 12a ) 

u'l = B(u' ± ) = u± ha - 2, (5.12b) 

4ac = «± = A(U±) = - U± Fba + ^ (5-12c) 

where, as indicated by the notation, it^ = Ao B o A(up ) are the phase-space parameters of the first epoch of the 
following era (which oscillates between the a and c walls, starting on a). Note that the transformations appearing in 
the era composition Ao B o A do not include the effect of the first wall b in (|5.6[) . Indeed, the collision on b would 
be (conventionally) included in the composition of transformations appearing in the previous era (which oscillated 
between c and b, ending on b). By contrast, wc conventionally include the effect of the collision on the last wall of 
an era in the composition of transformations associated to this era (e.g. the last transformation in (|5.12[) represents 
the last collision, on the wall a of the era (|5.6[) ) . The number of transformations which are composed during an era 
Ek(x, y) is equal to the length k of the era, i.e. to the number of arrows (or of epochs) in the diagram (|5.6|) of the era. 

The successive transformations, 

u -> —u -> u - 2 -> —u + 2 -> u - 4 -> ... (5.13) 

that appear in a long era of big-billiard oscillations between the a and b walls differ from the standard BKL result for 
oscillations between a and 6, namely: 

u^u-l^w-2^w-3^u-4— s-... (5.14) 

However, this difference is only due to the fact that the big-billiard description is keeping track of an information 
that BKL did not wish to keep track of: namely, the precise order between the Kasner exponents pi,p m (in the 
notation of BKL) associated to the 'oscillating' diagonal metric components a 2 ~ t 2pi , b 2 ~ t 2pm . As an unordered set 
{pi,Pm} = {pi,P2} (with pi < p2 < P3 as in section HTT1 above) . Indeed, we see in the list of 'Kasner transformations', 
Table HI labeling the possible permutations of the Kasner exponents, that the transformation 

u' = k 5 (u) = -u- 1, (5.15) 



4 As wc shall discuss below, the initial value of the m+ alone suffices to determine the explicit form of T among A, B, C. 

5 The child game called 'hopscotch' in English is called 'marelle' in French and 'campana' in Italian. 
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B a b 


-00 < u + < -1 


< u~ < oo 


Bba 


< U + < OO 


— oo < u~ < —1 


Bac 


-1 < u + < 


< < oo 


Bbc 


-1 < u + < 


— oo < u~ < —1 


Bca 


< u+ < oo 


-1 < u~ < 


B c b 


— OO < u + < — 1 


-1 < u~ < 



TABLE II. Hopscotch Court 



which corresponds to the permutation between p\ and P2, maps each one of the apparently discrepant values of u 
(namely —u, —u + 2, etc) in the big-billiard sequence (|5.13[) into the corresponding usual BKL one (|5.14[) . To wit 

fcg(— u)=u— 1; fcs(— u + 2) = u — 3; etc. (5.16) 

Let us now clarify what is the shape of the 'hopscotch court', i.e. the part of the (u~ , u + ) plane which parametrizes 
the dynamics of the big billiard. This full hopscotch court is naturally divided into six separate 'boxes':(l) a box, 
say B a b, parametrizing the epochs going from a to b, (2) a box, say B bai parametrizing the epochs going form b to a, 
etc, when considering the other ordered pairs (x, y) with (x, y) £ {a, b, c}. The precise boundaries of the box B xy are 
easily obtained by requiring that, in the Poincare half plane, there exists an oriented half circle (orthogonal to the 
boundary) crossing the walls x and y in that order. For instance, it is easily seen that the box B ab is defined by the 
inequalities 

B ab : < vT < +oo, -oo < u + < -1. (5.17) 

The inequalities defining all the boxes B xy are gathered in Table llll and the corresponding regions in the (u~,u + ) 
plane are represented in FigEJ Two important remarks concerning these boxes are: (i) all the boxes have a rectangular 
shape, and (ii) the union of all the boxes (together with their boundaries) does not cover the full (it~,w + ) plane. 
More precisely, the domain of the (u + , u~) plane which does not parametrize any epoch is the union of the following 
three 'vacuum boxes' 

V a : < u~ < +oo, < u + < +oo, (5.18a) 
V b : -oo < ir < -1, -oo < u + < -1, (5.18b) 
V c : -1< u~ < 0, -1< u + < 0. (5.18c) 

For instance, we illustrate in Fig. [S]a long era of epochs oscillating between b and a (see the points marked 1, 2, 3, 4, 5) 
which starts (point 1) in the upper right part of the box B ba , and then jumps successively from B ba toward B a b and 
back until it 'exits' by terminating in B bc (point 5). Then the (u - , u + ) point will jump from B bc to £? c t>, as part of a 
next era of the Ek (be) type. 



We have seen above that the (bijective) applications A, B, C (Eqs. (|5.2[) . (|5.3[) . (|5.4p ) corresponding to the collisions 

on the walls a, b, c, respectively, leave invariant the two-form Co^ed' ^' c l- l|4.11|) . that we shall simply denote in the 
following as 

„ du + A du~ ,„ „. 

UJ = 2 TT ^2- ( 5 - 19 ) 

(?v — u y 

Therefore the 2-form uj defines an oriented measure on the full hopscotch court that is invariant under the hopscotch 
discrete map T defined above. We would seem to be able to straightforwardly apply the tools and results of ergodic 
theory to our full hopscotch game. In particular, we know that this hopscotch game must be ergodic in the sense 
that it cannot leave invariant a subdomain, having a non-zero measure with respect to the form u, Eq. (|5.19p . of the 
full hopscotch court. Indeed, our hopscotch game is a projection of the billiard motion within an ideal triangle on 
the hyperbolic plane %2- If the projected billiard motion could leave 'unvisited', for an infinite 'time', a continuous 
subdomain of the full hopscotch court, this would be inconsistent with the fact that the billiard dynamics on an ideal 
triangle on H2 is known (since the classic work of Hcdlung and Hopf) to be ergodic, and therefore to visit the full 
3-dimcnsional phase-space (u~,u + ,s) which lies 'above' our hopscotch court. 

However, there is a catch in that, contrary to what is usually assumed in most investigations of ergodic theory, the 
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full invariant measure of our (projected) phase-space, i.e. the integral of the two-form lu, Eq. (|5.19j) . on the full 
hopscotch court defined in Table [TT1 is infinite). Indeed, the form (|5 . 1 9[) is singular along the line u + = u~, as well as 
at infinity ~ \u~\ ~ oo, where the integral of to is logarithmically divergent. We see on Fig S3 that the singular 
line u + = u~ lies mostly in the excluded ('vacuum') part of the hopscotch court. However, this line touches the 
boundaries of the court around the points (u~,u + ) = (0,0) and (u~,u + ) = (—1,-1). In addition, the boxes B a b 
and Bba extend at infinity, where J lu diverges logarithmically. It is also easily seen that the integral of lu over B ca 
and B ac diverges logarithmically near (u + ,u~) = (0,0), and that the same is true for J lu near (u + ,u~) = (— 1, — 1). 
The three points (0, 0), (—1, —1), (oo, oo) correspond to the three 'cusps' of the ideal triangle on H2- Indeed, another 

(2) 

way to understand why the integral of lu is infinite is to remember that Lu red = to can also be written (in 'Birkhoff 
coordinates') as dl A d(sma) = cosadl A da, where / measures the length of the boundary of the billiard, and where 
a, — 7r/2 < a < 7r/2, is the angle between the normal to the boundary and the velocity vector. The integral of sin a 
yields a factor 2, while the integral over dl yields the total length of the boundary of the billiard. In the case of the 
ideal triangle, the length diverges logarithmically at each corner. 



Before discussing other issues concerning the invariant measure lu, Eq. (|5.19[) . in the 2-plane (u~,u + ), let us note 
that, by marginalizing the variable u~ , we can deduce from to an invariant measure for the dynamics of u + alone. We 
already noticed that the hopscotch map T acts diagonally on (it - , u + ) (i.e. separately, and actually in the same way, 
on u~ and u + ). We warned the reader above that, in spite of this diagonal action, one needs to keep track of the 
action of T on the two variables (u~ , u + ) in order to determine the full, two-sided sequence of collisions corresponding 
to the billiard dynamics taking place within the chamber of the big billiard. However, if one ignores the variable u~ , 
and only considers the action of T on u + , it is easily seen that the sole knowledge of the initial value of u + suffices to 
determine the explicit expression of T (among A, B or C) and therefore all the future values of u + . Indeed, a look at 
Fig. [5] shows that there are three, and only three, possible cases: (i) if u + belongs to the interval [— oo, —1], the next 
collision will be on the b wall so that the action of T is u' + = T(u + ) = B(u + ); (ii) if u + G [—1,0], the next collision 
is on the c wall, so that T{u + ) = C(u + ); and (iii) if u + E [0, +oo], one has T(u + ) = A(u + ). Therefore, once we know 
u + , we can uniquely determine all its T iterates. [Reciprocally, it is easy to see that the knowledge of the initial value 
of u~ suffices to determine all the past values of u~, i.e. all its T _1 iterates.] In other words, if we simply ignore the 
variable u~ , the map T defines a dynamics for u + alone, which is an unquotiented version of the usual BKL dynamics 
on the single variable u recalled in the Introduction. [Remember that the BKL variable u actually coincides with our 
variable u + .] These remarks show that the unquotiented generalization of the BKL w-map defined by u' + = T(u + ) 
will admit an invariant one- dimensional measure w(u + )du + obtained by marginalizing (i.e. integrating upon) the 
variable u~ in the two-dimensional measure (|5.19[) . Explicitly, we can then define w(u + )du + as 

^du^lJ^du^J^J^, (5.20) 

so that 

w(u + ) = ^e b T * +v (5-21) 

Here, w^(it + ) denote the various boundaries of the integration domain on the u~ axis, and tb the associated signs. 
As one sees on Fig. [6j these boundaries are piecewise-constant functions of u + . For instance, when u + > 0, there 
are two boundaries: u~ iin = — oo (with e m j n = —1), and M„ lax = (with e max = +1). This leads to w(u + ) = \/u + . 
When —1 < u + < 0, one has four different boundaries (on each side of the vacuum domain in the middle of Fig. [S]). 
Finally, we can conclude that the unquotiented, big billiard dynamics of the variable u + considered separately leaves 
invariant the measure w(u + )du + where 



if < u+ < +oo, w(u + ) = — , (5.22a) 

u + 

if -l< U +<0, w ( u +) = ^—-^ = - 1 (5.22b) 

u + + 1 u + u + (u + + 1) 

if - oo < u + < -1, w(u + ) = — — . (5.22c) 

The existence of this invariant one-dimensional measure for the (unquotiented) BKL dynamics is not well known, and 
the explicit expression of w(u + )du + has, as far as we know, never been written down before. 

It was indicated in Fig. [5] that u + £ R should really be considered as a coordinate on the Kasner circle, i.e. the 
manifold of solutions pi,P2,P3 of the two Kasner constraints. Note that one could also parametrize the Kasner circle 
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FIG. 6. Billiard phase space in the u u parametrization: the epoch hopscotch court. The B xy regions are sketched, and filled 
with different colors (shades of gray) . 



by an angle 9 (with the usual 2-7T period). In addition, one can require that 9 be equal, say, to at the ab corner 
(u + = +oo), to 2n/3 at the be one (u + = — 1) and to 47r/3 at the ca one (u + = 0). With these requirements, and 
remembering that the disk model of Fig. [2] is related with the half-plane model of Fig. [3] (in which u + appears 
as a natural coordinate along the absolute) by a fractional linear transformation between the complex coordinates 
z = u + iv of Fig. 02 and £ = x + iy of Fig. [21 of the form 

C = 2z+1 ~^ , (5.23) 

the transformation between u + and 9 (as defined above) is given by (with £ = ie t6 on the Kasner circle, i.e. the 
absolute) 

- w - U+ + 5 ~ l ^ (5.24) 



2 ' 1 2 



The expression of the invariant measure w(u + )du + , Eqs. ()5.22j) . in terms of the angle 9 on the Kasner circle (in the 
disk representation) is then easily obtained by using Eq. (|5.24p . It is found to only depend on the angular distances 
between 9 and the two big-billiard "corners" surrounding it. Let us denote the angular location of the "corner a&", 
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between the walls a and b, sec Fig. [2] by 9 a b, and similarly for the angular location 9b c for the corner be, and 9 ca for 
the corner ca. These three angles correspond to u + = +00, —1 and 0, respectively, and with our chosen normalization, 
take the values 9 a b = 0, 9b c — 27r/3 and 9 ca = 4n/3 or — 2n/3. With this notation, the invariant measure w(6)dd on 
the Kasner circle is given, when 9 a b < 9 < 9b c (i-e. when 9 is on the "negative" , or "shadow" , side of the b wall) by 

w{9)d9 = ^ , e _ e f , _ e (9 ab <9<9 bc ). (5.25) 
4 sm % ab sm bc 2 

Its expression in the two other intervals 9bc < 9 < 9 ca and 9 ca < 9 < 9 a b is obtained by cyclic permutations 
abc — > feca — > ca&. Note that this measure is invariant around the middle points, e.g. {9 a b + 9bc)/2 for the interval 
[9 a b,8bc\ of Eq. (|5.25|) . and is logarithmically divergent as 9 tends to the extremities. For instance, as 9 — > 9 a b = 0, 
J e w(9')d9' ~ J„ d6' /9 1 ~ In 1/9 in keeping with the logarithmic divergence ~ ln(— u + ) of the corresponding u + interval 
(|5.22b ) as u + —> —00, given the link 9 ~ — v / 3/u + deduced from Eq. (|5.24[) as u + —> ±00. 

Let us note finally that it is straightforward to check directly the invariance of the measure w(u + )du + defined above 
under the explicit transformation laws A, B, C defined in Eqs. (|5.2[) , (|5.3[) . (|5.4p . However, to do that one must note 
that, contrary to the two-dimensional map T acting on the (u~ , u + ) plane (which is one-to-one), the one-dimensional 
action of T on the real line of the variable u + is no longer one-to-one, but rather two-to-one. [Indeed, the preimagc 
of a certain u + located within one of the three intervals [—00, —1], [—1,0], [0,+oo], can lie in either one of the two 
other intervals.] In such a case, one must remember [30j that the definition of the invariance of a measure fi under 
T is that, for any measurable set U, n{U) = /.t(T _1 ?7), where T -1 is the pre-image (rather then the image) of U. 
When applied to an infinitesimal interval / = u + ± ^du + , one must then take into account that T -1 / consists of two 
separate infinitesimal intervals. 

The invariant measure w{u + )du + (or w(9)d9) on the Kasner circle is not normalizablc. Indeed, we pointed out that 
its integral diverges logarithmically near each one of the three corners of the big billiard (i.e. near u + — 00, —1 and 0, 
or 9 = 0, 2tt/3 or Att/3). Therefore, the invariant measure of both the two-dimensional map T and its one-dimensional 
restriction are not normalizablc. It can seem strange to have an invariant measure of a projected phase-space which is 
infinitely large, while the invariant measure of the original, unprojected, phase-space was finite. Indeed, the invariant 
Liouville measure of the billiard on an ideal triangle, namely 

Q (3) = duAdv^dp 
v 

where the angle /3 (0 < /3 < 27r) parametrizes the angular direction of the unit velocity vector, integrates to the 
product of the finite area of the idea triangle J J dudv/v 2 = 2ir and of J d/3 = 2ir. As proven in full generality (for any 
time- independent Hamiltonian system) in the Appendix, see Eq. (|A. 13|) . the energy-shell-reduced Liouville measure 

is simply equal to the reduced product of the symplectic measure f^j^ oc (^rcd) A< '" _1 '' by ds, where s denotes a 
phase-space coordinate which is canonically conjugate to the Hamiltonian (so that ds/dt = 1 along the Hamiltonian 
flow). In the present case, where n = d — 1 = 2, this yields 

nf ] =uj{u-,u+) Ads (5.27) 

where s measures (when considering a unit velocity billiard) the hyperbolic length along the billiard orbit. This agrees 
(modulo an unimportant factor) with the result (25) of [23|. In terms of the phase-space coordinates (U,V,9), or 
equivalently (u + , u~ , 9^ of Eqs. (|4.7j) - (|4.9j) . we have ds = d9j sin#, so that 

, 1 N 9 _ 1 11 — cos 9 . . 

s[u ,u ,9) = so(u , ir) + In tan - = sq(u ,u H — hi -. (5.28) 

2 2 1 + cos 9 

Here, we can choose the u ± -dependent integration constant as we wish. For instance, we can choose it so that, for 
any given (u + ,u~), s(u~ ,u + ,9) varies between and some maximum value, say a(u~,u + ) as 9 varies between the 
starting wall and the ending wall along the oriented geodesic defined by (u~ , u + ). With this choice, we see that the 
full three-dimensional big-billiard phase-space has the shape, in (u~ , u + , s) coordinates, of a 'slab', above the (u~ , u + ) 
hopscotch court, of varying thickness < s < a(u~ ,u + ). Using Eqs. (|5.28[) and (|4.7[) . one can express the thickness 



6 The unit (Euclidean) velocity vector (vp u ,vp n ) = (cos /3, sin /3) = (sign(V) sin 9, cos 0) so that we have a link of the form /3 + 7r/2 = 
— sign(V)6 which ensures that the inequality < < 7r, together with the fact that V can be either positive or negative, corresponds to 
an angular direction f) of the velocity vector varying over a 2n range. 
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a(u , u + ) of this slab in terms of the values of the u coordinates (in the Poincare model) of the starting and ending 
walls along the geodesic (say u stai t{w ,u + ), u en ^(u~ , it + )), namely 



In 



^cnd 



In 



u start — u 



II 



u start 



1 , ("end —U )(u+ - M start ) 

- in - r — 

2 (u + — Ur, 



id)(«start - u ) 



(5.29) 



It is then easily checked that a(u~ , u + ) tends to zero near each corner of the billiard, thereby ensuring the convergence 
of the Liouville measure (|5.27p . i.e. the convergence of J J oj(u~ ,u + )a(u~ ,u + ). [For instance, near the corner 
u~ — > oo, u + — > oo] (with Ugtart = 0, M on d = — 1 or the reverse), the result (|5 . 29[) yields a thickness a(u~,u + ) ~ 

\ ("end - Ustart)(w+ 1 - UZ 1 ) -> 0. 



The ergodic theory of transformations preserving a measure on an infinite-measure space (or infinite ergodic theory) 
is an active field of current mathematical research in which, however, there are many less concrete general results 
than for the case of finite measure. In order to be able to avail ourselves of the usual ergodic theorems (such as the 
equality between the 'time average' and the 'measure average'), it is useful to transform the problem onto another 
one exhibiting a finite measure. Several different strategies are possible for doing so. 



As a first strategy, we could lift the full big billiard hopscotch game back to the hyperbolic billiard it came from. 
This would mean considering the ergodic properties of functions on the three-dimensional phase-space (it + , u~, s), 
with the finite Liouville measure (|5.27j) . In this case, we would be considering a continuous Hamiltonian flow, so that 
the relevant ergodic theorem would assert that, for almost every phase-space point x, 

lim I f dT>f{F T ,{x)) = (5.30) 

T^oo T Jo J fl{x') 

Here f(x) is a (measurable) function on phase-space, Ft(x) denotes the Hamiltonian flow over a timxjfl T, and n 
denotes the relevant, finite measure, i.e. fi = Eq. flOTl) . Then, if we were interested in the ergodic properties 
of phase-space functions f(x) = f(u~,u + ,s) that do not depend on s, i.e. on functions f(u~,u + ) that live on the 
(u + ,u~) hopscotch court, we can conclude from (|5.30|) (using the fact that ds/dT = 1 along each geodesic segment) 
that 

Um En=av(T n (u-,u+))f(T n (u-,u+)) = j^{u-,u+)a{u-^)f{u-^) 
n^L E^(MT™ («-,«+)) fu(u-,u+)o-(u-,u+) 

where o~(u~ ,u + ) is the thickness (in the s direction) of the phase-space slab 'above' the point (u~,« + ), and where 
T denotes the billiard map, i.e. the discrete hopscotch map transforming any (u~,u + ) parametrizing one Kasner 
epoch, into the values (u'_,u' + ) parametrizing the next Kasner epoch. As explained above, T is equal to A, B or 
C depending on the wall on which the considered geodesic segment will collide. The notation T n denotes the n-th 
iteration T° To ... oT, i.e. a composed transformation of the type of Eq. (|5.9|) . Note how the continuous time 
average of Eq. (|5.30[) has reduced itself (for functions depending only on u~ and u + ) to a discrete-time average, i.e. 
to a discrete sum (|5.3ip involving the successive iterates of the hopscotch map. However, the continuous-time origin 
of (|5.31[) is recalled through the occurrence of the 'weights' a (T n (u~ ,u + )) involving the successive thicknesses of 
the phase-space slabs, encountered along the billiard trajectory in (u~,u + ,s) space. [In (u~,ii + ,s) space the billiard 
motion becomes a so-called special flow [30l |. i.e. a combination of uniform motion in the 'vertical' s direction, with 
< s < a{u~ ,u + ), with discrete jumps in (u~,u + ) and in s (back to zero when it reaches cr(u~ ,u + )).] 



A second strategy for reducing the problem to a discrete map having a finite measure is to follow BKL in lumping 
together the epochs into eras, and to focus on the statistical properties of eras rather than epochs. In order to do 
this, we need to know on which subregions of the full hopscotch court, Fig. [51 each type of era E*(x,y) must start, 
[here, as above, (x, y) denotes an oriented pair of walls, and E„ denotes the union of all E^s, i.e. an era of arbitrary 
length k = 1, 2, 3, ... ] This is straightforwardly obtained by using the transformation rules A, B, C discussed above. 
For instance, the sub-region say F a b (where F stands for 'First') of the hopscotch court corresponding to the start of 
an era of the E* (a, b) type must come from a c wall, and include the process c — > a — > b. Using either some simple 



7 See, e.g., |31| for an entry into the literature on infinite ergodic theory. 

8 As recalled in Section III Al the appropriate time variable for the hyperbolic billiard is the T-time of Eq. 112.251 1. 
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TABLE III. Starting and ending subregions 



geometric reasoning, or working with the algebraic relations defining the transformations A, B and C, Eq. (|5.12l) . 
one finds that the F ab subregion is the rectangular subdomain of the B a b box defined by the inequalities 

F ab : < u~ < 1, -oo < u+ < -1. (5.32) 

The full set of inequalities defining the six possible starting subregions F xy , with x,y <S {a, 6, c} are listed in Table 
IIIII For compactness, we also indicate the six possible subregions on which an era of the type E*(x, y) can end. They 
are denoted by L xy (where L stands for 'Last'). If we consider the overlap domain F xy f] L xy between the start and 
the end of some (x, y)-type era, it must correspond to an era E\{x, y) of length k = 1, i.e. containing only one epoch. 
For instance, we see on Table iLLTl that the intersection Ff, a H-k&a corresponds to the small box 

F£ a : — 2 < u~ < — 1, < u+ < 1. (5.33) 

The box F^ a . Eq. (|5.33[) , describes the starting domain, in the (u + ,u~) plane, of all the one-epoch eras of the 
&a-type. More generally, it is not difficult to write down the inequalities defining the starting domains, say , of all 
the /c-epoch eras (with k = 1, 2, 3, ...) of the (starting) xy-type. They are given by intersecting the full F xy with the 
condition 

n xy {u+) = k (5.34) 

where n xy (u + ) is the integer- valued function listed in the second column of Table ITVl which yields the length of 
the era starting at some given point (u~,u + ) in phase-space. This leads to the 'era hopscotch court' of Fig. [7] 
which represents the six era-starting domains, and their division in A:-cpoch subregions F£ y . Note that the function 
n xy (u + ) giving the length of each era depends only on u + , and not on u~ . This corresponds to the fact that on 
Fig. [7] all the boundaries between the F£ boxes are horizontal. We have also indicated in Fig. [7J the special points 

(u-,u+) = (-0-1,0), (u~,u + ) = (-0,0+1), (u~,u+) = (-1/(2 + 0), -1/(1 - 0) (where 0= (y/E - l)/2 ~ 0.618 
denotes the 'small' golden ratio) corresponding to the simplest periodic hopscotch orbit corresponding to the infinite 
succession of one-epoch eras (|5.7[) . [There exists also the 'time-reverse' version of (|5.7p . namely a — > b — > c — > a... 
which jumps between F^ b — > F b \ — > F} a — > ....] 

For completeness, we have also indicated in Table HVl the discrete sequence of values of it^j, where 1 < m < n xy 
describing the successive epochs "contained" within an era that starts from some it* 6 F™ y (one example of such 
sequence of epochs was drawn in Fig. [B]). In Table HVl u m7 . y denotes either u^^^ or u0 xj/ [we indeed recall that the 
discrete hopscotch map T acts on a diagonal manner on u + and u~: u' + = T(u + ) and u'_ = T(u~)]. Moreover, within 
some era E n ^ v (x, y) starting with an epoch of the x — > y type, roughly half of the epochs contained in E n *y {x, y) are 
of the x — > y type (namely those corresponding to m = 1,3,5,...) while the other half are of the y — > x type (those 
corresponding to m = 2,4, ...). 



9 The notation [x] for x g IR denotes the usual integer part of x when x > (e.g. [n] = 3), and — [— x] < when x < (so that [— n] = —3). 
We did not find useful to introduce other definitions of the integer part (e.g. the 'floor', 'ceiling', or 'Hurwitz' ones). 
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TABLE IV. Epoch Hopscotch 



B. Era hopscotch dynamics 

If, following the spirit of Belinski, Khalatnikhov and Lifshitz, we focus on the discrete dynamics of successive eras, 
we can consider a hopscotch game based on the 'era hopscotch court' represented in Fig. i.e. the six era-starting 
domains F xy (further divided in sub-boxes labelling the length of the era) . The resulting discrete era-transition maps 
Tern (mapping the (w~~,u + ) point of the first epoch of an era to that of the first epoch of the next era) will be 
obtained by composing the individual epoch-transition maps contained in the considered era. E.g. an era E^{a,b), 
i.e. a — > b — > a — > b, would correspond to T era = B o Ao B. Depending on the parity of the number n xy of epochs 
contained in the considered era E n *y(x,y), the era-transition map T e va will map the initial starting rectangle F xy to 
a next starting rectangle, say F' x , y , where the labels x' and y' are fully determined by the knowledge of (x, y) and of 
the parity of n xy (i.e. whether it is even or odd). The explicit rules giving F' x , y , for each F xy are given in the first 
columns of Table [Vj In addition, the explicit form of the corresponding era-transition map, 

Up, = Tn"v [Up J , (5.35) 

transforming the phase-space point u ^ of the first epoch in some era E n x V (x, y) into the phase-space point u» 

of the first epoch in the next era E nX i v ' (x' ,y') are explicitly given in the last column of Table IVl For instance, if 
we consider up ba = (u F ,<j>), <f> ~ 0.618 denoting the small golden ratio as above, we shall have (from Table HV]) 
n ba = [<f>] + 1 = 1, which is odd, so that (from Table IVf the next era will be F' ac , and the new starting phase-space 
point in F' ac will have as coordinates (from the last column of Table IIV[) 

u F ^ = -Up, u F ^ = -cj>. (5.36) 

We have proven above that the 2-form u>, (|5 . 19[) was left invariant by each individual wall collision transformation A, 
B or C (and more generally by any symmetry transformation of H^)- Therefore, L) will be invariant under the era- 
transition maps 7~ n *y , which are certain compositions of n xy wall-collision transformations, e.g. 7~E 3 (a,b) — B o Ao B. 

A crucial property of the era-Hopscotch court, Fig. [7J is that the integral of uj over the era court is finite. Indeed, the 
points (0, 0), (— 1, —1) and (oo, oo) leading to the logarithmic divergence of the w-measurc of the full epoch hopscotch 
court of Fig. [5] are well separated from the six F xy era rectangles of the era hopscotch court. For instance, the region 
at infinity of the Fb a rectangle is —2 < u~ < — 1, u + — > +oo, which leads to convergence for Jw(u _ ,h + ). 

We are therefore in the usual conditions for applying the results of ergodic theory on a finite-measure space. In 
other words, after many iterations of the era map 7^ r a we can consider that the phase-space point (u~,u + ) behaves 
in a stochastic manner, described by a 'probability measure' equal to u>/I, where I is the integral of uj over the era 
hopscotch court of Fig. [JJ Note that, in the present context of the iteration of a discrete map T, which is ergodic 
and admits an invariant measure uj, the 'probability' measure is uj itself, and its meaning is that the ratio J_^uj/ ^uj 
(where A C T> is a subregion of the full domain T> of the era hopscotch court) yields the n — > oo limit of the fraction 
n^/n of the number of eras n_4 spent in A among the n first iterates of an arbitrary initial phase-space point. In other 
words, the word 'probability' refers here to a limiting era-frequency. We shall explicitly compute some probabilities 
in the era-hopscotch dynamics below. 
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FIG. 7. The era hopscotch court in the (u~,u + ) parametrization. The six starting boxes F xy are sketched, and filled with 
different colors (shades of gray), according to Fig. [6] as explained in Table Hill For each starting box F xy , the starting sub-boxes 
corresponding to fc-epoch eras (as described in Table HV)l are indicated for k = 1,2,3; namely F xy , F xy and F xy . The three 
(u~ ,u + ) points of the simplest periodic orbit are denoted as asterisks. Please note that, for typographical reasons, the -F^'s 
are in this figure indicated as Fkxy. 



VI. SYMMETRY-QUOTIENTING THE BIG BILLIARD 



As already mentioned, the basic a, 6, c system, Eq. underlying the big billiard dynamics is formally invariant 

under the six-fold group of permutations of the 3 letters a, b, c, say S3. This group S3 is the symmetry group of 
the (ideal) triangle, in the Lobachcvsky plane, represented on Fig[TJ It comprises (when seen in the disk model) two 
rotations by ±27r/3 (that exchange the corners among themselves), and three reflections with respect to the lines 
bisecting the corners (that permute two sides among themselves). Note also that the action of the six elements of S 3 
on the boundary of the disk correspond to the 5 Kasner transformations given in Table (JTJ) , together with identity 
transformation, say kr, (with v! = kr,(u) = u). 

We can use the symmetry group S3 to quotient the dynamics of the big billiard. There are (at least) two ways of 
thinking about this quotienting. One way would be to consider a kaleidoscopic version of the big billiard dynamics in 
which the single "moving ball" of the billiard, shown e.g. in Fig. @J is augmented by its 5 (generically distinct) images 
under S3 . This leads to a billiard game where 6 (symmetry related) balls simultaneously move within the same billiard 
table, and (simultaneously) bounce on its bounding walls. The phase-space points of this kaleidoscopic billiard is a 
symmetry orbit of the original (single ball) phase-space point, i.e. an unordered set of (at most) six (two-dimensional) 
g l 's and six (two-dimensional) p^s (restricted by the condition g^PiPj = 1). [Some of these variables are allowed to 
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TABLE V. Era Hopscotch 



coincide when the original ball crosses one, or several, of the fixed sets of 53 (bisecting lines).]. 
A second way to look at the quotiented dynamics is to replace the latter kaleidoscopic phase-space point 
{9(1) > 9(2) > •"' 9(6)'^i > — >Pi ) by its unique representative, say <j* cp , within a fundamental domain of 53, together 

with its corresponding momenta p* op . For instance, one can use as fundamental domain the small billiard, with sides 
G,B,R in Fig. [5] [We shall use this relation when considering the small billiard dynamics below.] When passing 
from the continuous billiard dynamics to the discrete billiard map from an epoch to the next epoch, the quotienting 
of the big billiard leads to a quotiented version of the epoch hopscotch game of Fig. |6] For instance, the kaleidoscopic 
version of Fig. [6] would replace each (u~,u + ) point indicated there by six symmetry-related points, i.e. u and its 
five transforms under the Kasner transformations of Table U namely /c2(m ± ), fc3(it ± ), fe^ir^), fcs(u ± ). These 

six points would then simultaneously "jump" , after a (six-fold) collision on a gravitational wall, to their next six-fold 
positions in the (u~,u + ) plane. For instance, the S3 orbit {l 7 , l 11 , l 111 , 1 IV , l v , 1 VI } of the point 1 in Bh a (Fig. H]) 
would jump onto the new 53 orbit {2 Z , 2 7/ , 2 /7/ , 2 /l/ , 2^, 2 177 }. In the alternative, fundamental-domain, version of 
the quotiented dynamics we could replace each S3 orbit in the (w - ,u + ) plane by its unique representative located 
within, say, the box B^ a . [Indeed, the six boxes B xy of Fig. [5] are exchanged under S3.] In that view, the discrete 
quotiented big billiard would become a map from the box B ba onto itself. For instance, in the example shown in 
Fig. IH1 the initial point 1 would first jump to the point 2' (midway between the points 1 and 3), then to a point 4' 
(midway between 2 and 5, and belonging to Bb a ). The next epoch would be the image of the point 5 £ Bh c which 
belongs to Bb a ■ As suggested by our description of the example of Fig. El one finds that each era gets quotiented into 
a succession of (u~ , u + ) representative points within B\> a which lie on a straight (Euclidean) line of slope +1 (as the 
line passing through the points 1 and 3 in Fig. [5]). More precisely, if we denote by (up, B^ a the first epoch of a 
quotiented era, the i?b a -representative of the version of the considered era is made of the points 

(u F ,up'); (u F -l,u^-l); (u F - 2, up - 2); (u F - [up], up" - [up]) (6.1) 

obtained by successively subtracting 1 both from u F and ut As indicated, the length of the era is simply given by 
k = [lip] + 1, so that the last point of an era is reached when its u + coordinate is betweenO and 1: u\ = u\ — \op\ = 
{u F }. Note, however, that the vT coordinate of the last epoch is given by u£ = u F — k + 1 = u F — [u F ], so that it 
depends both on u F and on u J (while the sequence of the u + values depends only on the starting value of of u + ) . 
Then, it is easily seen that the next epoch (i.e. the first point of the next era) will be (when mapped back to Bt> a by 
S3) 

( 1 + -1, + * 
\u F -[up] up -[up] 

If we ignore the u~ coordinate, we see that the law giving the successive values of the u + coordinate coincides with the 
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law found long ago by BKL (when using the BKL„ >0 definition of an era, as discussed above) namely: Up = x + k—1, 
Up — 1 = x + k — 2, down to = x, with k = [u~p] + 1 denoting the length of the era. This shows that the BKL 
discrete dynamics of the variable u is obtained by: (i) quotienting our more complete hopscotch dynamics by the 
permutation group S3, and (ii) ignoring the u~ coordinate and identifying the BKL variable u with u + . Note again 
that this link between the hopscotch dynamics and the BKL dynamics is particularly simple if one uses the BKL u> o 
definition of an era, rather than the BKL u> i. 

The S^-quotiented hopscotch dynamics, i.e. the discrete dynamics mapping Bb a onto itself, defined by Eqs. (|6.f [) and 
(|6.2p . constitutes a two- variable generalization of the BKL map, say 

(u~,u + ) 6 B ba ->• (u~ ,u + ) = % a (u~,u + ) e B ba . (6.3) 

Like the full (unquoticnted) epoch hopscotch dynamics, the quotientcd discrete map, T ba , (|6.3p leaves invariant the 
measure oj However, like in the unquoticnted case, the integral of to on the domain (and image) Bba of the map Tba 
is infinite. 



Before discussing the obtention of a finite- measure discrete map associated with % a , let us note that, as we did 
for the unquotiented hopscotch dynamics, we can also consider the action of the quotiented billiard map Tba on the 
single (BKL- like) variable u + . It is defined as 

u+ -> u+ - 1 -> ... -> u+ - [u+] -> 1 - 1 -> ... (6.4) 

U F ~ \ U f\ 

In other words, it is just the usual BKL map on u recalled in the Introduction (in its BKL u >o version). Note the 
appearance of —I in the definition of the new "first u + " of the next era. This takes into account that the unordered 
triplet of Kasncr exponents of {uj} < I is identical to that of {up} -1 and should therefore not appear twice in the 
sequence of transforms of u + . 

As in the unordered case discussed above, we can obtain an invariant measure for the one-dimensional version, say 
T b a °f the quotiented billiard map by marginalizing over iT the two-dimensional invariant measure Wba (Le. the 
restriction of ui to the region Bba)- We use the same definition as above, Eq. (|5.20[) . The difference is that, now, 
u + is restricted to be in the interval < u + < +00, and the boundaries of integration over u~ are u~ lin = —00 and 
u mas = — I (where u~ ax differs from its previous value). This leads to the invariant one-dimensional measure 

du + 

Wba(u + )du + = — — — . (6.5) 
vp~ + 1 

It is again easy to check directly that Wb a (u + )du + is invariant (in the sense recalled in Section lTvl and in the Appendix) 
under the discrete map T b ^ [one needs to take into account that TSP is two-to-onc so that the preimage of an 
infinitesimal interval u + ± ^du + consists of two infinitesimal intervals.] Note that if, instead of considering the 

definition above of T b ^\ one considers the "standard" BKL u> i map (where u remains always > I and decreases by 
units of 1 until the value I + {u}, after which it jumps to l/{u}), the expression of the invariant measure reads 

du 

WBKL u>1 (u)du = — (6.6) 

As far as we know, the results in Eqs. (|6.5p and (|6.6p have not been explicitly discussed before in the literature. 
It should be noted that the one-dimensional measure Wb a (u + )du + differs from the restriction to Bba of its unquotiented 
analog discussed in Section [V] above. [Indeed, we had before w(u + )du = (du + /u + when u + > 0.] This is an effect of the 
quotienting which means that one must fold back onto Bb a the symmetric images that were elsewhere (and notably 
in B ca ). Finally, we note that, as before, the invariant measure w{u + )du + is not normalizable. Indeed, it diverges 
logarithmically when u + —> +00. Note, however, that it converges at the lower boundary u + — > 0. [The same is true 
of the invariant measure u>bkl u>1 (u)du = 1/u of the standard BKL u-map, with I < u < +00]. 



Let us now come back to discussing the two-dimensional quotiented map Tba, acting on Bb a - To get a discrete map 
preserving a finite measure, we need to consider the quotiented analog of the era hopscotch dynamics. In fact, we 
obtain an era hopscotch dynamics simply by ignoring the intermediate epochs and focusing on the map transforming 
the (quotiented) first epoch of an era (up,Up) into the first epoch (u F ,Up) of the next era. We shall denote this 
quotiented era map as T. When using, as we did above, a representative of the S3 orbit within the B ba box, the 
quotiented era map T is a map of Fb a onto itself. We recall that Fb a is the domain of the first points of eras starting 



29 



A-, 


ba 
U K = 


—U a b — 1 


A-2 




— (Uac + 1)/Uac 


ka 


ba 
U K = 


—Ubc/{Ubc + 1) 


A-i 


ii ba — 


l/Uca 


A'4 


„.ba 

U K — 


-l/(uc + 1) 



TABLE VI. The Kasner maps for each starting box of the big billiard table 



as b — > a. It is the rectangular domain — 2 < u F < —1, < Up < +oo (see Fig. [7]). The explicit expression of the 
map T is given by Eq. (|6.2I) . i.e. 



Ju f = +— — - !■ (6-7) 



1 

up - [u FS 

This map leaves invariant the restriction of the two-form u> to the domain i 7 ^, i.e 



u F = 2 ?i A % (6.8) 

( U F - u f) 

By contrast with the original measure uj on the full hopscotch court, this restricted measure has now a finite integral, 
namely 



cjF = 21n2 (6.9) 

We shall also refer to the map T as being the Chernoff- Barrow- Lifshitz- Khalatnikov- Sinai- Khanin- Shchur map 
for the big billiard, or CB-LKSKS map in short. This one-to-one map between two variables was introduced in [9I) 10 I. 
while Ref. and 11 1 showed how such a two- variable map appears as a completion of the original BKL analysis, 



when keeping not only the original BKL variable u + (parametrizing pi, pi and ^3), but also a variable 8 related to 
the amplitude of oscillations of the a, 6, c metric variables during an era (the precise relation between our variables 
and those used in [l(| and [ll| will be given below). The CB- LKSKS map is one-to-one over its domain Fh a . Note, 
however, that the original map T on the era hopscotch court, whose quoticnting leads to T, is such that both the 
image, TF^a, and the pre-image T _1 i 7 b a , of F^ a are the union of an infinite number of rectangular domains belonging 
to two different F xy boxes. For instance, the image 7~F), a is of the form Pb c U P QC , where Pb c is the union of an infinite 
number of disconnected rectangles contained within i 7 ^, while P ac is the union of an infinite numbers of disconnected 
rectangles contained within F ac . 

Note that up~ and Up play asymmetric roles in the T map. Indeed, for all the starting boxes F xy , the boundaries of 
the subdomains corresponding to a given era- length are horizontal (see Fig. [7]). Therefore, whether a given rectangle 
in Fb a is T-mappcd into one or several rectangles only depends on the range of up~ independently of the range of Up. 
As a consequence, the 'mixing' character of the T map is essentially contained in the u + direction. Actually, we see 
by differentiating Eq. (|6.7[) that 

As —2 < Up < — 1 and [vip] G N, the denominator (u F — [uj]) 2 is always strictly larger than one. Therefore, 
|dTup| /dup < 1, i.e. the T map is contracting in the Up direction. 

Finally, let us clarify the link between the CB-LKSKS map as defined above and the unit-square map given in [Tol | 
and [ll|. The range of the variables u~ , u + in the above-defined CB-LKSKS map is the infinite vertical rectangle Fb a , 
i.e. —2 < u~ < — 1, < u + < +00. By contrast, the statistical analysis developed in [l(| and [ll[ is based on two 



The definition of the map in [j| and our definition match when one relates the variable denoted as x in [j| , which we shall call Xqb , 
with our Up by XnB = — l/( n F + !)• For the sake of completeness, let us also note that the names of the Kasner exponents pi and P2 
are exchanged in 
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variables x + , x_ , defined in the unit square, i.e. < x+ < 1 , < x + < 1 . 

Let us consider the following transformation mapping Fb a into the unit square < x + < 1 and < .t_ < 1: 

x+ = (6.11a) 

Up + 1 

= — up - 1. (6.11b) 
In terms of the unit-square variables (x-,x+), Eqs. ()6 . 11 1) . the CB-LKSKS map reads 

Tx + = {±}, (6.12a) 
1 



T. 



X- 





1 


X- + 


x +. 



(6.12b) 



consistently with the results of [TT| - As the unit-square transformation in (|6 . 1 1 1) does not consist of applying the 
same fractional linear transformation to both variables u + and u~ , the two-form uj is not invariant under (|6.1ip , but 
becomes 

dx+ A dx- 

(1 + X+X-Y 

The choice of Bb a (and Fb a ) as representative boxes for the quoticnted big billiard dynamics docs not affect the unit- 
square results (|6.12p and (|6.13[) . Had we considered another representative box B xy , and defined a correspondingly 
modified version of the unit-square transformation (|6.1ip . we would have ended up with the same results (|6.12p and 

(ESI). 



VII. SOME PROPERTIES OF THE SYMMETRY-QUOTIENTED DYNAMICS 

In this section, we shall start from the -Ff, a -box version of the CB-LKSKS map (or, simply, T map), Eq. (|6.7[) . and 
study some of its properties, recalling, when needed, some results obtained by previous authors and focusing on new 
results. 



A. Probabilities 



Let us indicate how, in our set-up, one can compute the probability for an era to contain a given number of epochs. 
The probability P ni for an era to contain a number n\ of epochs, 1 < n\ < +00, is proportional to the integral of 
the two- form cop (|6.8p over the relevant region of the u F ,Up space. In the considered case of a given number, m, of 
epochs, the relevant region is simply the box , as defined above, i.e. the domain —2 < Up < —1, m — 1 < up~ < n\. 
To normalize this probability, we must then divide by the integral of cup over the full domain Fba, a s given in Eq. 
(|6.9p . 2 In 2. As a result, we obtain that the probability P ni for an era to contain a number ri\ of epochs as 

1 n2j ni _ 1 7_ 2 ( u f- u f) ln2 rii(n 1 +2) 

This agrees with the result of BKL (which was obtained from the stationary probability distribution (|1.4p of the 
Gauss iteration map). 

For instance, the probabilities for the length of an era to take the values n\ = 1, 2, 3, 4, 5 are 

Pi = 0.4150 (7.2a) 

P 2 = 0.1699 (7.2b) 

P 3 = 0.0931 (7.2c) 

Pi = 0.0589 (7.2d) 

P 5 = 0.0406, (7.2c) 
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whose sum is 

Efe=i p k = 0.7775. Hence, 77.75% of the eras have lengths smaller or equal to 5, and actually 58.49% 
of the eras have lengths smaller or equal to 2. This shows that most eras have rather small lengths. The asymptotic 
behaviour of P ni , as m — s- +oo is 

P ni * ^ (7.3) 

with c = In 2 ~ 0.69315. Therefore, the probability to have n\ > N\ is asymptotically given (for large A^i) by 
P{n\ > Ni) ~ which decreases rather slowly as N\ increases. In other words, though most eras have rather 
small lengths, from time to time eras with unbounded lengths can arise. The rather slow decrease of P ni as ni 
increases implies, in particular, that the mean value of n\ is infinite (being given by the logarithmically divergent 
series Sni>i^"i n i)- Note, however, that a finite result is obtained if one considers the expectation value of the 
geometric mean of large sequences of independent era lengths, i.e. the exponential of the expectation value of lnni 
(known as the Khinchin number). This yields [IH 

exp p «i ln "i - 2.6854... (7.4) 

ni<l 

This result confirms that the "typical" length of an era is rather small. 

Let us now consider the computation of a probability of a more specific event (not explicitly considered by BKL): 
namely the probability P ni ,n 2 f° r an era of length 1 < n% < oo to be followed by an era of length 1 < n 2 < oo. This 
probability is obtained by integrating the form (|5.19[) over the appropriate range of the variables Up,u F . This range 
is determined by Table pvp (as in the previous case), and by the properties of the CB-LKSKS map: combining the 
two information, the range n\ — 1 + n ^ +1 < u^ < ri\ — 1 + is obtained. The range of Up is determined as in the 
previous case. As a result, we obtain 



-1 r n 1 -l+- 1 



p 1 / A,r I 1 ~' n2 dU F - 1 In f ("i"2 + l)(nirc 2 +ni+n 2 +2) \ (7 K\ 

*ni,n 2 — In 2 / ua F 1 / + „ -x 2 ~ In 2 111 V (fiin 2 +ni + l)(nin 2 +n 2 + l) I ' ^'-°> 



2 J ni _i+_± FT (Up- Up) 



Note that this probability is symmetric in n\ and ri2- The previous probability P ni is (as it should) recovered by 
summing P Tll .„ 2 over all the values of 712. 

n 2 — 00 

P = -±- V P =^-\n( ( " 1+1)2 "| (7 6) 

n 2 = l 

Note that in both cases we were considering events that depend only on Up so that the result involved marginalizing 
the variable Up, i.e. integrating the (normalized) oj form over the complete range of Up. And, indeed, integrating the 
(normalized) two-form ujf over the complete range of Up (—2 < Up < —1) yields the one-form 

1 f^ 1 du F du+ _ 1 dup ^ ^ 



ln2j_ 2 {u+ -Up) 2 ln2( u + + l)( u + + 2) 

In2 T+x 



which yields the Gauss distribution (jl.4p . i.e. w(x)dx = when parametrizing up~ g]0,+oo[ by x g]0, 1[, such 



that u F = i — 1 (with < x < 1). 
Note that the (integrable) invariant one-dimensional measure (|7.7p differs from the (non-integrable) invariant one- 
dimensional measure (|6.5p : while the former refers to the discrete dynamics of the u + value of the first epoch era, the 
latter refers to the discrete dynamics of the « + values of all the epochs. Note also that the probability P m to have 
an era of length n\ is obtained (in view of the link ri\ = [u~p] + 1) by integrating the measure (|7.7|) over the interval 
ni — 1 < Up < n\. 

It is interesting to note that the (pseudo-)random variables n\ and 712 (i.e. the lengths of two consecutive eras) 
are not independent of each other, because P ni n 2 7^ Pn 1 Pn 2 - However, the variables n\ and 712 are approximately 
independent statistical variables. Indeed, using the explicit expression (|7.ip and (|7.6p one finds that the ratio R ni n 2 = 
Pnmn /Pni Pn 2 j which would be, by definition, equal to one if n\ and 712 were independent random variables, takes values 
rather close to 1. For instance, R n ^ 0.8826, R 12 = R21 ^ 0.9985, R 13 = R 31 ~ 1.051 and R 22 ^ 1.008. Therefore 
the low values of n\ and n 2 (which are of greatest importance for many issues) are approximately independent. As 
concerns the large values of n\ and n 2l let us note that the asymptotic value of P ni n 2 is 
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where c = ln2 as above. This implies that the ratio R ni n 2 = Pnm 2 1 PniPn 2 ' 1S asymptotically constant, and equal to 
c m 0.69315. Having seen that the lengths of consecutive eras are approximately independent random variables, we 
expect that such an independence property will become more and more exact as one considers eras that are more and 
more separated. 

B. Continued fractions 

Let us briefly recall (from Q ) the usefulness of the continued- fraction representation of the variables u~ , u + in 
describing the effect of iterating the CB-LKSKS map. 
Any number < y < +oo, can be uniquely decomposed as 

y = [y} + M, (7.9) 

where [y] is its integer part, while {y} is its fractional part. This decomposition can be iterated by considering the 
decomposition (|7.9[) of l/{y}. This leads to the unique continued- fraction decomposition of y > 

M = «iH ; i — = [ni;n 2 ,ri3,...}, (7.10) 

where n\,n 2 ,ri3, ... are natural integers. In Eq. (|7.10[) we have introduced a notation for the continued- fraction 
expansion which distinguishes (by means of a semi-colon) the first integer m = [y\. In the case where n\ = (i.e. 
in the case < y < 1, i.e. y = {y}), we shall also use the notation {y} = [712,113,714, ...] (without semi-colon). The 
continued- fraction expansion contains a finite sequence of integers m, n 2 , ■■■ if y is rational, while it contains an infinite 
sequence of integers if y is irrational. 

With this notation, the continued-fraction expansionJ^l of the first (u~ ,u + ) values of a (quotiented) era starting 
(as above) in the box F ba can be written as 

Up = n ba - 1 + [n 2 , n 3 , n 4 , ...] = [n ba - 1; n 2 , n 3 , n 4 , ...], (7.11a) 
u F = -1 - [mi,m 2 ,m 3 ,...]. (7.11b) 

Here, we have denoted the first integer of the decomposition of as nb a — 1 = [up], 80 that n ba = [tit] + 1 denotes 
the length of the era starting with Up. In terms of these decompositions pertaining to the first era, we can write the 
first u~, u + values of the iV-th era (with N = 2, 3, ...), i.e. the (N — l)-th iteration of the T map 

1 u F ~ = n N - 1 + [n N+1 ,n N+2 ,n N+3 , ...], (7.12a) 

J N 1 u F = -1 - [njv,njv_i, ...,n 2 ,n ba ,m 1 ,m 2 , ...]. (7.12b) 

In other words, at each iteration of the CB-LKSKS map, the information about the length of the corresponding era 
is 'transferred' from the u + variable to the u~ variable. Note that the information contained in up' (i.e. the sequence 
of integers vn\,m 2 , ...) is progressively decaying in the "tail" of the iterates of Up. By contrast, the iterates of 
progressively uncover the "tail" of , thereby exhibiting the chaotic character of the Up dynamics which progressively 
amplifies smaller and smaller details of the continued- fraction expansion of the initial up'. 

C. Recovering information about the unquotiented dynamics 

Up to now, in this Section, we have discussed the quotiented dynamics (viewed within the representative box B ba , 
or F ba when considering the first epoch of an era). This quotienting has ignored the fuller information contained in 
the original, unquotiented era hopscotch dynamics, namely the precise corner, {xn, 2/jv} (with xn,Vn £ {a, b, c}) and 
orientation, (xn, J/^t), (i.e. xn ~> Un), of the first epoch of the 7V-th hopscotch era. For each starting box, say F xy , 
the specific Kasner transformation mapping the region F xy into the representative region F ba is given in Table I VII 



Here, we adopt the definition of integer part of a negative number given in Footnote [9] and we define the fractional part of a negative 
number accordingly, e.g. {— it} = —0.14.... 
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Therefore, if we start an unquotiented era hopscotch dynamics in some specific region F xy , we can first map it to the 
reference region Fb a by some specific Kasner transformation k xy to get its 6a-representative, up^ ha \ namely 

Then starting from [up ]& a we can define its continued- fraction decomposition, say 

uP^ ] = k xy u F ^ = n xy - 1 + [n 2 , n 3 , ...}. (7-14) 

Here, the so-defined integers n xv ,n2,n>3, ... give us the values of the successive lengths of all the eras that will evolve 
from the initial value ut . 

This reasoning shows how the knowledge of any era-starting values (up,Up) in the (u~,u + ) plane determine the 
subsequent era-length history. First, the location of up in the plane determines the initial corner xy. Second, 
the knowledge of the initial corner uniquely determines k xy (mapping it to Fba). And, third, the computation of 
k xy u~^ and of its continued- fraction expansion, determines all the era-lengths n xy , n%, TI3, .... However, we need 
more information if we wish to recover the full hopscotch dynamics from the simpler quoticnted hopscotch dynamics 
discussed above. Specifically, we need to recover information about: (i) what is the succession of corners (among the 
three corners of the big billiard) that will be visited, and (ii) what is the succession of the directions (clockwise or 
counter-clockwise) in which the oscillations within these corners will take place. Let us now show how one can recover 
this missing information from the knowledge of the era-starting values (up^ , ) in the (u~ , u + ) plane. 

Before doing so, let us establish some notation. Looking at Fig. [3J we shall say that an era is clockwise if its first 
epoch connects two billiard walls in the clockwise sense with respect to the unit disk, i.e. if it is either ba, ac or cb. 
In the other case (first epoch of the type ab, be or ca) we shall say that the era is counterclockwise. 
Now, we remark that the information about this direction of motion is contained in the determinant of the Kasner 
transformation k xy that maps each era-starting box F xy onto the reference box Fb a (as given in Table IVI| . More 
precisely, if the determinant of k xy , say D[k xy ] is equal to +1, the era F xy is clockwise, while if D[k xy ] is equal to 
— 1, the era is counter-clockwise. Given two era-starting regions F xy and F x > y i, it will sometimes be convenient to say 
that they are "parallel" if they have the same direction of motion (clockwise or counter-clockwise), and "antiparallel" 
in the other case. 



Le us now show how one can encode the information which is missing in the quotiented billiard in a pair p, r\, where 
p takes three different values, and rj two different ones. In more mathematical terms, p £ Z3 and r\ £ Z2, where Z n 
denotes the cyclic (multiplicative) group of order n. The values of p can be {1, exp(i27r/3), exp(— i27r/3)}, and are 
encoding rotations in the disk model of the billiard by the angles 0, 27r/3 or —2tt/3 respectively. The values of 77 are 
{+1,-1} and can encode the two possible "directions of motion" of an era (clockwise or counter-clockwise). 

Our aim is, starting from some initial (era-starting) position ( u Fj. y , u ~Fj. y )& Fxy, to determine the ordered corner 

xn,Vn within which the iV-th unquotiented T iterate, ~T N (up,Up) of (up,Up) will oscillate (i.e. the first epoch of 
the ( N + 1-th era) . We parametrize the ordered corner xm,Un by the pair pn,Vn where p^ is the rotation, in the disk 
model, mapping the initial (unordered) corner {x,y} into {xn,Un}, and where rjpf gives us the relative orientation 
between xn, Un and x, y. (i.e i]n = +1 if they are parallel, and tjn = —1 if they are antiparallel). 
One can iteratively build the values of p^ and t]n by using the following elementary facts: 

• an (intermediate) era F x i y i containing an odd (respectively, even) number of epochs is followed by an era F x ii y ii 
whose relative direction of motion is parallel (resp., antiparallel); 

• the rotation (in the disk model) between some intermediate era F x i y i and the following F x >> y n is equal to e +i27r / 3 
(resp. e - t2 ' n 7 3 ) if F x i y > contains an even (resp., odd) number of epochs. 

To exhibit the iterated effect of these elementary rules, it is convenient to define the following quantities (taking the 
values ±1): 

e j = (7.15) 

where nj denotes the number of epochs contained in the j-th era defined by the initial value of uj. With this notation, 
the Z3 x Z2 valued pair (pn,Vn) giving the rotation (with respect to the original era-starting domain F xy ) and the 
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relative "sense of motion" (parallel or antiparallcl to xy) corresponding to the ordered corner of T N (u F , Up) is defined 
by (with D[kp] = ±1 denoting as above the determinant of the Kasner transformation kp mapping (u F ,u F ) to Fb a ) 

PN = e ll>N , with e N = -D[k F } 2 §-(ei+e 1 e2 + ... + e 1 e 2 ...e N ), (7.16a) 
Vn = eie 2 ...ejv- (7.16b) 

Note that tjn can be rewritten as tjn = {—)® N , where 

N 

Q N = N + Y,n k . (7.17) 

fc=i 

Note also that the absolute sense of motion (clockwise or not) of the iV-th era is given by 

D[k F }r, N = £>[fc F ]e 1 e 2 ...e JV . (7.18) 



Given this result, we can now write, for any starting point (u F ,u F ), the explicit result of iterating N times the 
unquotiented T map, i.e. 

ut =T N u$ =T o ...oT oTu± , (7.19) 

in terms of the simpler action of the quotientcd T map, namely 

u F ^ N ^ N = kj yN T k xy Up xy . (7.20) 

[For brevity, we have denoted above the phase-space point (u F ,u F ) simply as u F ]. Here, k xy is, as above, the Kasner 
transformation mapping the initial era-starting box F xy to Ft, a: and k~}, yN is the inverse of the Kasner transformation 
mapping the ordered corner of T N up onto the "standard" ba corner. The transformation is determined from the 
above-computed values of (pN,T)isr). More precisely, the procedure determining k~^ yN is : (i) starting from k xy and 

the continued-fraction decomposition of k xy up~ = n\ — 1 + [712,71.3, ...], one determines the ej, Eq. (|7.15[) . and the pN 
and 77 at; (ii) then p^ determines the rotation between xy and the final corner xpj,yN, and tjn determines whether 
this corner is "parallel" or "antiparallel" to the initial xy; (iii) finally, knowing the ordered corner (x N , y N ), Table IVll 
determines the transformation k x N v n that maps it onto the ba corner. 



Let us note the dissymmetric roles of u + and u~ . In the above construction, it was the knowledge of the initial 
value of ?i + which allowed one to recover the full information about the future evolution of the unquotiented dynamics. 
The situation would be different if we wanted to describe the past unquotiented dynamics. In that case, it would be 
the continued-fraction of the ba transformation of the initial u F that would encode the needed information. 



To make the above construction more concrete, let us end this subsection by working out an explicit example. 
We consider (for simplicity) an initial era-starting box of the ba type. For the convenience of the reader, we list in 
Table IVlTl the concrete meaning, for this case, of the six different values of (pn, Vn) in determining the ordered corner 
of the iV-th iterate of the initial point. Let us for instance consider tj^ + 1 = [m; tt-2, •■•] of the form u F + 1 = \/2, 

i.e. u Fb + 1 = [1; 2, 2, ...], and consider the second iterate T 2 u Fb . In that case, xy = ba and k xy = fco (the identity), 
and therefore D[k xy ] = +1. As n\ is odd and 712 even, we easily find that pi — 1 and 772 = — 1. This shows that the 
second iterate (i.e. the third era, if we count the initial one) is of the ab type. The explicit expression of T^Up is 
then 

T*u% a =k~jT\% a , (7.21) 

where k a b = &5 in the list IVII We have thereby reduced the computation of the iteration of T to the simpler 
computation of the iteration of its quotientcd version T. 



D. On the anisotropic behavior of the unquotiented big billiard 



The aim of this subsection is to highlight one interesting feature of the unquotiented big billiard that is lost in its 
quotientcd description: its anisotropy, i.e. the fact that, after each given era (taking place in some corner, with some 
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r?M 


M M 

x y 


1 


1 


ba 


1 


-1 


ab 




1 


ac 




-1 


ca 


■ 4 


1 


cb 


■ 4 


-1 


be 



TABLE VII. The exit possibilities for a sequence of eras starting with Fb a - The rotation pu is considered modulo 2n. 



sense of motion for the first epoch) the next era has more probability to take place in a specific ordered corner, namely 
a corner obtained from the previous one by rotating it in the same direction as the first epoch in the disk model, and 
keeping the same sense of motion for the first epoch (by ±27r/3). For instance, if the first era is, say, of the ba type, 
the following era has more probability to be of the ac type rather then the be one (which is the other possibility) . 

Indeed, the general formulas (|7.16l) above show that, when N = 1, i.e. after one iteration, the second era is obtained 
from the first by applying the rotation p\ = e 01 , with 6\ — ~^-D[Kp\ei, and that its relative sense of motion is 
?/i = e\. Here, t\ = (— ) 1+ni is determined by the parity of the length ri\ of the first era. On the other hand, the sign 
ei, which determines both p\ and 771, is a statistical variable whose probability distribution is determined by that of 
ni, i.e. by P ni , Eq. ()7.1[) . Among the two possible values of ex, the most probable is ei = +1, corresponding to ri\ 
being odd. Indeed, this probability is obtained by summing (|7.1[) over n\ = 2k + 1, k € N, and reads 

P(C1 . +1) . M _ £ .J. h ( (2 ^+ 2 f +3) ) = ^ * 0,5,5 

The complementary probability that t\ = — 1 i.c that n\ be even, is 

In other words, we have a strongly anisotropic behavior after one era: if the first era is, say, of the ba type, the 
following era will be of the ac type in 61.15% of cases, and of the be one in only 34.85% of cases. 

Let us now see what happens after two iterations. Fixing for simplicity the initial ordered corner to be ba, the 
ordered corner after two iterations can be of four different types: 

• if (61,62) = (+1, +1) (corresponding to 62 = — 4tt/3, 772 = +1 ), it will be cb; 

• if (ei, 62) = (+1, —1) (corresponding to 62 = 0, 772 = —1 ), it will be ab; 

• if (£1,62) = (— 1, +1) (corresponding to 62 = 47r/3, r/2 = —1 ), it will be ca; 

• if (ei, £2) = (— 1, —1) (corresponding to 62 = 0, 772 = +1 ), it will be ba. 

The probabilities corresponding to each one of these cases is easily computed from the probability distribu- 
tion (|7.5p . for ni,7i2 (remembering that ei = (— )" 1+1 , £2 = (— )™ 2+1 ). For instance, P{e\ = +1,^2 = +1) 
= J2k 2 >OM>o P2k ^+ 1 > 2k ^+ 1 = ^odd.odd is found to be 



Similarly, 



P{ei = +1, e 2 = +1) = Podd.odd =s 0.4199. (7.24) 



P(d - +1, e 2 = -1) = Podd, cvc „ = J2 p 2fc 1+ i,2fc 2 ^ 0.2316, (7.25a) 

fe 2 >l,fci>0 

P(C1 - -1, £2 = +1) = Pcvcn.odd = J2 P ^,2k 2 + 1 = Podd,cvc„ 0.2316, (7.25b) 

fc 2 >0,fci>l 

P(ei = -1, £ 2 = -1) = Povcn.cvcn = XI P^fei ,2fc 2 ^ 0.1169 (7.25c) 

fe 2 >l,fel>l 
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Again we see a strong 'anisotropy' among the various possibilities. In particular, the most probable case is again that 
corresponding to applying the two rotations (by ±27r/3) in the direction indicated by the first epoch, keeping the same 
sense of motion. The calculation gets more involved for higher iterations. One expects that, after many iterations, 
the memory of the initial ordered corner will get lost, and that one will end up with asymptotically equal probabilities 
in any one of the six possible ordered corners. Note, however, that the anisotropic behavior we are discussing here 
will continue to be present locally: after each era (whether or not it is much "later" than the initial era), the next era 
will take place in a preferred corner with respect to the previous one. 



Let us briefly discuss periodic orbits in the big billiard with a focus on the differences between the notion of periodic 
orbit in the unquotiented billiard, and the corresponding notions either in the quotiented billiard or in the BKL map 
(acting solely on u + ). 

Given any discrete map, say T, acting on some space X, a periodic orbit is a set of successive T images of a point 
x e X, say {x,Tx, . . . ,T m ~ 1 x}, such that T m x = x. The (minimal possible) integer m is called the period of the 
discrete map T. As the CB-LKSKS T map is a quotiented version of the full hopscotch map T, it is easily seen that 
any n-periodic orbit of T will automatically "descend" to a corresponding periodic orbit of T. However, the period 
of the corresponding T orbit might be a divisor of m. On the other hand, it is a priori possible that periodic orbits 
of T could not be "lifted" to periodic orbits of T. To study these two issues (the change in period from T to T, and 
the possibility of lifting periodic orbits from T to T) let us start from some given n-periodic orbit of T. 

A first issue that should be discussed is the relation between periodic orbits of the BKL map Tbkl (i-e. the restriction 
of the two-dimensional map T to the one-dimensional map u'< = Tbkl(m+)), and periodic orbits of T in the (u~ , u + ) 
plane. The periodicity condition of Tbkl involves only one condition, namely u + = Tg^ L (w + ) = T'"(u+), while the 
periodicity condition of T looks much more restrictive as it involves two separate conditions, namely u + = T (u+) 
and U- — T (u-, (we recall that the action of T on u+ only depends on w+, while its action on u_ depends, 
in addition, on the integer part of u+). However, we have seen above that the T map was always contracting in 
the u~ direction (see Eq. (|6.10j) V Therefore, we expect that the iterated effect on any starting value of u~ of the 
finite collection of maps indexed by the various values of [u+] in the periodic orbit T(u~ , [u + ]) will converge to some 
corresponding fixed orbit of u~ values. 

Let us then start by an arbitrary m-periodic orbit of the one-dimensional BKL map: T^ Kl ^{u + ) = u + . It is 
well-known, [33j . pij . and evident from the explicit form of the action of T on the continued- fraction expansion 
of u + , that any such m-periodic orbit is parametrized by the special values of u + that admit a (regular) periodic 
continued-fraction expansion (cfe) of the type 



By well-known theorems going back to Euler and Lagrange all such values of u + + 1 are quadratic irrational numbers, 
i.e. irrational real roots of quadratic equations of the form ax 2 + bx + c = with integer coefficients and a positive 
discriminant b 2 — 4ac (that is not a perfect square). The simplest example of such pcriodic-cfe numbers is the (large) 
golden ratio u+ + 1 = $ = [1;1, 1,1, 1,. . .] = (a/5 + l)/2 ~ 1.618 (so that w+ = $- l = 0= (a/5 - l)/2 is equal to 
the small golden ratio <fi = $ — 1 ~ 0.618). The second simplest examples are of the type u + + 1 = [n; n, n, n, . . .] = 
{s/n 2 + 4 + n)/2, with some integer n > 1 and are sometimes called "silver ratios". Note that in the cases of the 
golden ratio, or of the silver ratios, the period m of the BKL orbit is m = 1. 

Let us first show that any m-periodic orbit of the one-dimensional BKL map Tbkl gives rise to a unique correspond- 
ing m-periodic orbit of the two-dimensional T map. This follows from the explicit expression (|7.12[) of the iterated 
action of T on an arbitrary starting value of u~ , written as Eq. (|2.32[) . By repeatedly iterating the (w + -dependcnt) 
action of T on u~ one sees that the information contained in the initial value of u~ is lost in the receding tail of the 
cfe of I u so that the sequence of values of u tends to a fixed orbit that is entirely defined by the periodic cfe of 
u + . More precisely, the limiting value of u~ , which pairs with the given u + to define a two-dimensional m-periodic 
orbit of T is given by 



VIII. PERIODIC ORBITS 



u + + 1 = [ni;n2,ri3, . . . ,n m ,ni,n2,ns, ■ ■ ■] ■ 



(8.1) 



(i + O = [ 



(8.2) 
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Having shown that any m-periodic orbit of the one-dimensional BKL map T, u + — > Tbkl(w + ), i.e. any (regular) 
periodic cfe of the type, Eq. (|8.ip . uniquely determines a corresponding m-periodic orbit of the two-dimensional T 
map, (u + , u~) —> T(u + , u~), we now discuss the issue whether any m-periodic orbit of T can be lifted to some periodic 
orbit of the full, unquotiented billiard. This question can be answered positively by studying the m-iterated action of 
the unquotiented map T on the initial point (it~ , u + ) € Bi, a of a m-periodic orbit of the quotiented map T. The fact 
that T (it~,u + ) = (u~,u + ) means that the ^-symmetry orbit of T m (u~ ,u + ) coincides with that of (u~,u + ). This 
means that there exists a particular Kasner transformation fc* (which depends on m and on the considered periodic 
orbit of T) such that 

T m (u-,u+) = k,(u-,u+). (8.3) 



The set of six Kasner transformations is a realization of the S3 permutation group (of order 3! = 6). In fact, this 
permutation group consists of the identity, 3 transpositions [(12), (23) and (31)], and 2 cyclic transformations [(213) 
and (321)]. We recall that the order of a particular group element, such as fc*, is the smallest integer p such that 
fc* = k . As a transposition is of order 2, and a cyclic permutation, (123) or (321), of order 3, we see that the order 
p of fc* must be equal to p = 1, 2 or 3. Therefore, by iterating (|8.3[) . we get 

T m P(u-,u+) = K{u-,u+) = (u- u+), (8.4) 

and mp will be the smallest such integer. In other words, (u~,u + ) is the initial point of a periodic orbit under the 
unquotiented billiard map T, with period pm, where p — 1, 2, 3 is the order of k*. 

Actually, the specific value of p can be algorithmically derived from the value of the Z2 x Z% pair {rim , p m ) corresponding 
to T m . More precisely: (i) if i] m = —1, then p = 2, independently of the value of p m ; (ii) if r\ m = 1 and p m = 1, then 
p = 1; while (hi) if r/ m = 1 and p m =/= 1, then p = 3. (Note that, when m = 1 =>■ p rn =/= 1.) 



We have therefore proven that any periodic orbit of the quotiented map T (or, even, any periodic orbit of the 
one-dimensional BKL map Tbkl) can be lifted to a periodic orbit of the unquotiented big billiard map T. Note 
that this property extends to the corresponding continuous billiard motion in the unquotiented big billiard (simply 
by considering the geodesic segments corresponding to all the (u~,w + )'s belonging the periodic orbit). To make this 
general result more concrete, let us consider a particular example. The simplest periodic orbit of the quotiented 
billiard is that given by the golden ratio, namely 

u + + l= [1; 1,1,1,.. 
-(«- + l) = [1,1,1, 

"Downstairs" its period is m = 1. However, the T transform of the above "golden-ratio" point (u~ , u + ) = ( — 1 — </>, (j>) 
is given by reflection in the a wall, i.e. by the matrix A of Eq. (|5.14p . so that 

Tu + = -u + = -<(>, (8.6a) 
Tu- = -vT = ! + (/). (8.6b) 



V5+ 1 
2 

V5- 



5.5a) 



(8.5b) 



The transformation v! = —u is not one of the Kasner transformation. However, in keeping with the general reasoning 
above one can use the hxed-point property of this golden-ratio periodic orbit (namely 0(1 + 4>) = 1) to rewrite the 
r.h.s.'s of the above equations as 

Tu+ = --i- = -^^r = k 4 (u+) 7 (8.7a) 
1 + u + + 1 

TW = - = — — — = ki{u~) (8.7b) 
<p u + 1 

where k±{u) = — l/(u+l) is one of the Kasner transformations of Table|U Therefore, the specific Kasner transformation 
fc* corresponding to the particular "golden-ratio" periodic orbit of T is fc* = k±. The latter Kasner transformation 
correspond to the element of S3 realizing the cyclic permutation (j>i,P2,P3) — > (j>2,P3,Pi) of Kasner exponents. The 
order of such a cyclic transformation is p = 3 : k 4 o k 4 o fc 4 = identity. This shows that the golden-ratio initial 
conditions above define a periodic orbit of the unquotiented billiard of order pm = 3x1 = 3. We recover the periodic 
orbit of Eq. (|5.7[) made of three successive one-epoch eras between the middle of the three successive gravitational 
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walls a, c and 6. 

Note that, in the general case of a starting value for u + of the type (|8.ip . an m-periodic orbit of the quotiented 
billiard will contain m eras containing, successively, nx, Ti2, ■ • ■ , n m epochs (so that it contains nx + n.2 + • ■ • + n m 
epochs in all), while the lift of this periodic orbit onto the unquotiented big billiard will contain pm eras, containing 
p(nx + n<x + . . . + n m ) epochs in all. 

Finally, let us note that the billiard periodic orbits discussed here are the projection down to hyperbolic space H2 of 
Lorentzian-billiard motions in /3-space which arc not periodic there. Indeed, the spatial metric gij(T, x) corresponding 
to these dynamics is expressed in terms of the /3's, rather than the projected 7's, say (for the diagonal Bianchi IX 
case of relevance to the big billiard) 



with 



and 



^(T,x) = ^e- 2 ^( T ) e ^(x) ej a (x), 

a 

p a (T) = p(T) 1 a (T) 7 
p(T) ~ exp(cT) , (with c> 0) 



(8.8) 

(8.9) 
(8.10) 



where we used the result [17[ , [2l| that A = In p is (asymptotically) a linear function of the coordinate time T defined 
in Eq. (|2.25p . A periodic orbit of the big billiard is such that 7 a (T + nP) = 7 a (T) for n € N and some period V 
in T-timc. This periodicity "downstairs" in T-time does not correspond to a periodicity of the metric coefficients 
gij(T, x). It does not even correspond, as one might have thought, to a discrete self-similar symmetry of the metric 
(i.e. gij{T + nP,x) = A n <?y (T, x)) but to a rather different discrete transformation under which the "scale factors" 
a a (T) = e~P ( T ) (i.e. the BKL variables a,b,c) change as 

a a {T + nV) = [a a (T)] x " (8.11) 

with 

A = e cV > 1 . (8.12) 

As the 7°'s are confined by the big billiard walls to remain (non strictly) positive, the (periodic) scale factors a a (T) 
stay < 1, and we see on Eq. (|8.11[) that (apart when they collide on a gravitational wall, where the corresponding 
scale factor becomes equal to 1) the scale factors tend to zero super exponentially with the number n of periods. 

IX. SMALL BILLIARD 

Up to now we have been discussing the big billiard, with three walls a, 6, c making up an ideal triangle in hyperbolic 
space (see Fig. [5] or Fig. [3]). This billiard corresponds to the dynamics of the diagonal Bianchi IX model, i.e. the 
a, b, c system of BKL. However, as recalled above, in the most general (non-diagonal, inhomogeneous) case, the use 
of an Iwasawa decomposition of the spatial metric, as in Eq. (|2.4j) . leads to a closely related but slightly different 
billiard, namely the small billiard made of one gravitational wall, and two symmetry walls. As in Fig. [5] or Fig. [3J we 
shall denote the (partial) gravitational wall as G (for gravity or green), and its symmetry walls as B (for blue) and 
R (for red). More precisely, in the notation of the Poincare model of Fig. [3l 

• the G wall is the portion u = 0, v > 1 of the gravitational wall a = 0; 

• the B wall is the portion u = —1/2, v > V3/2 of the symmetry wall a = b: 

• the R wall is the portion u 2 + v 2 = l,— l/2<u<0of the symmetry wall a = c. 

Our aim in this section is to relate the dynamics within the small billiard to the dynamics within the big billiard 
studied above. Somewhat surprisingly, though the two billiard tables are closely related, the two corresponding 
dynamics cannot be straightforwardly mapped among themselves. 
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BG 


U < — 1, U + > U a 

-1 < u~ < -1/2, u+ > up 


BR 


U~ < — 1, Up < U + < U a 


RG 


— 1 < U~ < —1/2, U a < U + < Up 
-1/2 < 7i~ < 0, ?i + > U a 


RB 


-1/2 < u" < 0, m+ < 
< u~ < 1, u a < it+ < 


GR 


> 1, U a < M + < «g 


GB 


< U _ < 1, U + < It a 
> 1, U + < U,£j 



TABLE VIII. The regions of the u + u plane, where the dynamics of the small billiard takes place. 



A. Dynamics of the unquotiented small billiard 



We start by considering the dynamics within the three walls G,B,R of the small billiard, without introducing any 
extra quotienting. Indeed, the small billiard table is already a fundamental domain of the six-fold symmetry group 
S3 of the big billiard acting on I-L2, therefore one could a priori expect that the small-billiard dynamics be equivalent 
to the quotient of the big-billiard one by S3 (as studied above). Actually, this is not the case. We shall find that 
the small billiard dynamics is not equivalent to the quotiented big billiard dynamics. The basic reason for this non 
equivalence is that the small billiard tabic is obtained by quotienting only the configuration space (q space) of the big 
billiard dynamics by S3, while the S3-quoticnting we considered above was done in phase-space {q,p space). When 
considering on its own the small billiard dynamics it is natural to introduce analogs of the notions introduced (by 
BKL) in the big billiard context. First, we shall define an epoch of the small billiard as a geodesic segment (i.e. a 
Kasner motion) connecting two successive walls. For example, a B — > G epoch (or, for short, a BG epoch) is an epoch 
starting from the blue (B) wall and ending on the green or gravitational (G) wall, etc. 



The dynamics of successive epochs of the small billiard is similar to that of the big billiard. For instance, the 
following sequence of epochs 

R^G^B^G^R (9.1) 

corresponds to a succession of "collisions" on the G, B and G walls for a dynamics which started on the red wall and 
returned on it. As in the big billiard case, it is convenient to parametrize each epoch by a point in the (u~ , u + ) plane, 
with it + , resp. u~ , parametrizing the end, resp. beginning, of the extended geodesic corresponding to the considered 
epoch. A first difference with the big billiard case is that the regions of the u~u + plane which describe the small 
billiard dynamics are quite dissimilar to the corresponding big-billiard regions drawn in Fig. 6. 

They are drawn in Fig. [SJ they comprise six allowed regions, labelled as BG, BR, RG, RB, GR and GB, and a 
large, connected forbidden ("vacuum") region which occupies the central part of Fig. 6 (between the upper hyperbola- 
like curves and the lower ones). The precise definition of the allowed regions is given in Table fVIIII where we used a 
short- hand notation for the following functions of vT : 



1 



11 



(9.2a) 



- s (9 - 2b) 

As we see the boundaries of the various allowed (and forbidden) regions are made of segments of hyperbolas: u + = 
(au~ + b)/(cu~~ + d). This is different from the big billiard case, where all the boundaries corresponded to horizontal 
or vertical lines (see Fig. [5]). In addition, the forbidden region of the big billiard (white domains in Fig. [5]) was 
made of three disconnected pieces. If we compare Fig. [5] and Fig. |S]we can roughly consider that the small-billiard 
u~u + picture, Fig. [SI is obtained by "morphing" the big-billiard one, Fig. 6, via a deformation where the B^a box 
becomes the BG one, B ca — > RG, Bb c — > BR, B ac — > GR, B c b — > RB and B a b — > GB, while the three disconnected 
forbidden regions of Fig. 6 "percolate" among themselves into the connected central forbidden domain of Fig. [5J This 
correspondence between the two pictures exists because, for instance, the extension to the big billiard of a small- 
billiard BG epoch corresponds to a ha epoch, etc. However, contrary to what one might have naively expected, it is 
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FIG. 8. The phase space of the small billiard in the parametrization. The regions -RG and RB are delimited by solid 

thick lines, the regions BR and GR are delimited by dashed thick lines. The Kasner intervals are indicated on both axes by 
thin gray lines. The subdomains corresponding to the first few eras of RG and RB are also sketched. 



not possible to find a globally-defined transformation u + = U(u~ , u + ), u~ = V(u~ , u + ) (leaving invariant the 2-form 
u, Eq. (|5.3[) ) that maps Fig. [5] into Fig. [S] The main obstacle to the existence of such a transformation is the fact 
that the forbidden regions of the small billiard comprise (say in the Poincare model) not only all the forbidden regions 
of the big billiard (geodesies in the Poincare half-plane that do not intersect the abc triangle), but, in addition, new 
forbidden regions (geodesies that intersect the abc triangle but not its GBR subtriangle) that were allowed before. 

Let us now briefly describe the dynamics of the small billiard, as seen in the u~u + plane, Fig. [5J It is similar to 
the "hopscotch game" associated to Fig. [5] Namely, each initial point (u~,u + ) in the occupied regions of Fig. [5] will 
"jump" to another position according to the following rules: 
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• if the point (u ,u + ) belongs cither to the BG region or the RG one (so that it corresponds to an epoch starting 
either on B or on R and ending by a collision on G) it will jump by a G-collision, i.e. by the transformation: 

G waU (u = 0) : m ± G^) = -u* (9.3) 

• if (m~, u + ) belongs either to the RB or GB regions, it jumps by the transformation: 

B wall (u = — -) : -> B{u ± ) = — - 1 (9.4) 

• if (u~ , u + ) belongs either to the BR or GR regions, it jumps by the transformation: 

R wall (u 2 + v 2 = 1) : u ± -> i?(u ± ) = l/u* . (9.5) 

Note that a point in some region will jump either to the YX or the FiT region (with {X, Y, Z} = {G, i?}). 
Note also that the "jumping rules" above are similar to, but different from, the corresponding A, B, C jumping rules 
of the big billiard, see Eqs. (|5.12|) . 

The small-billiard hopscotch game defined above leaves invariant the 2-form uj, Eq. (|4.11j) . However, as in the big 
billiard case the integral of w over the allowed regions of the small billiard is logarithmically infinite. As in the big 
billiard case, this logarithmic divergence comes (when writing w in terms of Birkhoff coordinates along the sides of 
the billiard table) from the infinite hyperbolic length of the G and B walls that meet on the absolute. This suggests a 
natural way of by-passing this divergence problem: to consider the (Poincare) return map of the small billiard on its 
red (R) wall, which is the only wall having a finite length. In other words, it is very natural, within the small billiard, 
to collect together all the epochs corresponding to bounces between the G and B walls into small-billiard eras, and 
to focus on the small-billiard era hopscotch dynamics which maps the beginning of such an era to the beginning of 
the next one. For instance, a small-billiard era comprising four epochs was indicated in Eq. (|9.1j) . For definitcness, 
we shall define the first epoch of a small-billiard era as the epoch which starts on R (for instance the leftmost R — > G 
epoch in Eq. (|9.ip ). Clearly, the era- hopscotch transformation Up = f(w^) mapping the coordinates of the first 
epoch in an era to the coordinates of the first epoch in the next era, being obtained by composing the individual 
G, B, R rules above, will be given by some diagonal fractional-linear transformation: 



v 



'f = T sb u F , u' F = T SB u' F , with T SB u = — ' " , (9.6) 

cu + d 

and will leave invariant the restriction to the R wall of the usual two- form u>, Eq. (|4.11[) . More precisely, the analog 
of what was before the "first-epoch" domains (such as the Fb a subregion of Bi, a represented in Fig. 7) become the 
"i?-lcaving" domains. Contrary to what happened in the big-billiard case, we do not need here to delineate these 
regions as subregions of the full small-billiard hopscotch court of Fig. [3] Indeed, by definition the era-starting region 
of the small billiard consists of the union of the RG and RB regions in Fig. [5] The RG region represents the 
first-epochs of the small-billiard eras that leave R towards the G wall, while the RB region represents the first-epochs 
of the small-billiard eras that leave R towards the B wall. They are indicated by solid thick lines in Fig. [H The 
era-hopscotch rule Tsb (|9.6p will map RGURB onto itself (mapping sometimes, say, RG to itself or to RB, etc.). Now, 
the small-billiard era-map Tsb leaves invariant w on a space (RG U RB) on which u has a finite integral. Therefore 
we are now in the good conditions for applying crgodic theory, and considering that to defines a probability measure 
on RG U RB. More precisely, one finds that 

l w = -lnf-V / -^ = -^2 (9.7) 
rg 2 2 \2J J RB 2 2 

so that the normalized probability measure on RG U RB is 

du+Adu- . , 1 1, /3\ 1, „ 1, n 

c— with - = -ln - +-In2 = -ln3. (9.8) 

(u+-u-) 2 c 2 V 2 / 2 2 

It would seem that, at this stage, we have obtained a small-billiard era-dynamics that is very similar to the 
quotientcd big-billiard era dynamics, i.e. the CB-LKSKS map T between the big-billiard first-epoch region Fb a to 
itself. An apparent qualitative difference is the fact that the small-billiard era-starting region RG U RB is made of two 
disconnected pieces. However, it is easily found that if we transform the RB region by the diagonal transformation 

u'± = -u± - 1 (9.9) 



au + b 
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(which leaves invariant uj), the RB region will be mapped onto a new region say RB' , of the u~u + plane which is 
contiguous to RG (along its left boundary). Therefore, we can simply replace the RG U RB domain by a thicker, 
connected domain RG U RB' (delimited in Fig. [8] by a dotted line) , and consider the dynamics of the (suitably 
transformed) era-hopscotch map Tsb from RG U RB' onto itself. At this stage, we have an era-hopscotch map which 
looks quite similar to the quotiented map T from F{, a onto itself. Moreover, if we look at the asymptotic region of 
long eras, i.e. the u + ^> 1 region of RG U RB' (where the RB' part came from the u + <C — 1 region of RB) we see 
that the shape of our new small-billiard first-epoch domain is simply 

-l<u~<0, u+>l (9.10) 

In other words, it is asymptotically rectangular and, modulo a simple shift of u~ by one unit, seems to coincide with 
the (exactly rectangular) quotiented big billiard domain Fb a - 



This asymptotic coincidence (modulo some suitable identifications) between the small-billiard red-return map and 
the quotiented big-billiard era- map T was physically expected because long eras of the big billiard (with u + 3> 1) 
correspond to many bounces between the b and a walls which clearly (see Fig. 5 together with Fig. 3) will correspond 
to roughly twice as many bounces between the B and G walls of the small billiard. Here, we have in mind folding 
some, say, a — > b — > a oscillation of the big billiard back onto the small billiard by introducing a B "mirror" in the 
middle of the ba corner, so that the big billiard bounces a — > (B) — > b — > (B) — > a (which do not "see" the B wall) 
become transformed in small-billiard oscillations of the type G^B^G^B^G, where we used the fact that, 
deep into the corner, a and G coincide. However, this asymptotic (large u + ) coincidence does not extend to the 
small-it + (short-era) region. In other words, there does not exist a w-preserving transformation of the u~u + plane 

that would map the RGURB >RGURB dynamics onto the Fb a > Fb a one. Indeed, if such a transformation 

existed the w-area of FR U RB (or equivalently RG U RB') would be equal to that w-area of Fb a . However, we have 

1 f 1 1 3 1 1 f 1 

-oj+ / -u = -ln- + -ln2= -In3^1n2 = / -u . (9.11) 
rg 2 J rb 2 2 2 2 2 J Fba 2 

However, this non equivalence does not mean that the small-billiard is unrelated to the big-billiard. [Actually, 
we shall discuss in the following subsection a well-defined, accurate relation between the two billiards.] It mainly 
means that the natural definition of eras in the small-billiard, as i?-return maps, cannot be identified with the natural 
definition of eras in the quotiented big billiard. It is true that, roughly speaking, the small-era definition (which 
means, when viewed in the big billiard, that an era ends when one crosses one of the three symmetry segments 
connecting the center of the disk to the middles of the three gravitational walls in Fig. [2]) signals the passage from an 
oscillation in one of the corners of the abc billiard to an oscillation in another corner. However, the problem is that 
the precise definition used in the big billiard of the transition between an afe-typc oscillation to, say, a foe-one cannot 
be equivalently characterized as the crossing of one of the three symmetry segments (similar to the R segment in 
Fig. 2). This non equivalence essentially stems from the fact that a big-billiard epoch (going to a gravitational wall to 
another gravitational wall) can, on its way, cross either 1, 2, or even 3 symmetry walls. And even if one considers the 
three symmetry segments of the type of the R one in Fig. 2, some big-billiard epochs can cross either 1 or 2 i?-type 
symmetry segments. This non constancy in the number of crossing of symmetry walls or segments prevents one from 
being able to define, in a uniform manner, a transformation between big-billiard eras, and small-billiard ones, which 
respects, say, the number of epochs during an era. 



In addition to this non equivalence in the definition of eras, a technical inconvenience of the unquotiented small- 
billiard era-hopscotch map defined above is that it is more difficult to find the generic, explicit expression of the 
small-billiard era map Tsb, Eq. (|9.6p . Indeed, while it was easy to define the small-billiard epoch hopscotch map, see 
Table HXl and Eqs. (|9.3[) , (|9.4|) . and (|9.5[) for any starting point u~u + , the corresponding definition of the small-billiard 
era map Tsb (as the return map on the R wall) has remained rather implicit. Given some (u~ , u + ) £ RG U RB one 
needs to iterate the epoch map a certain number of times (which depends on the starting point) to find the explicit 
expression of the red- return map Tsb- Actually, with some more effort it is possible to define Tsb nearly explicitly. 
Let us indicate how. First, we can delineate the small-billiard analogs of the F^ y boxes of the big billiard, i.e. the 
subrcgions of the RG and RB regions that will lead to (small-billiard) eras having some given number of epochs. 
Let us call RG n (resp. RB n ) the sub-region of RG (resp. RB) which provides a starting point for an n-epoch era. 
Here, as it is easily seen, n > 3 for RG and n > 2 for RB. [For instance, any i?G-starting era must go through 
R — > G — > B — » R before returning to R.} These regions are defined by the domains described in Table IIXI In this 
table the last column lists the equations of the four curves that delimit, for each n, the corresponding era-starting 
subdomains. These curves are given by equations of the type u + = f(u~), where the functions f(u~) are either the 
functions u a (u~), up(u~) defined in Eqs. (|9.2[) , or the following functions 
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RG", n odd 


U + > 1, M < — 1/2, (u a ,U,3, 


RG n , n even 


u+>i, u -<-i/2, (u,,,^;- 1 ,^ 1 ) 


RB n , n odd 


u+<0, t»->-l/2, (ucU^U^UZ) 


RB n , n even 


u+<o, u ->-i/2, {u a ,u p ,ur\u: +i ) 



TABLE IX. Small billiard hopscotch 



n _ 1 -2nu +2u +ra -2n+5 . 
b 2 -2«~+n-l ' 

, a 1 — 4re+7— 2nu~ -j-4u~ -j-?i 2 . 



-2-2«-+j! 



rrti 1 2nu 

u n — 2 



-2m +w--2rt+5 . 
2«-+?i-l ' 



f7' 



1 3+2nu_+2rr 

2 2u-+n 



Some of the corresponding regions, for n < 5, are represented in Fig. [5J As we see, contrary to the big-billiard 
case where the boundaries between the F£ y boxes where always horizontal (i.e. of the type u + = const), here the 
boundaries between RG n or RB n boxes are curved. This contributes to the difficulty of writing an explicit expression 
for the red-return map Tqb- Anyway, given these results, we can semi-explicitly define T$b by the following algorithm. 
Given some red-starting point (u~u + ), one must first find to which box, RG n or RB n it belongs by using Table llXl 
Then, knowing the type (RG or RB), and the length (n) of the era starting from (u~u + ), one can write T$b by 
composing n — 1 transformations of the G or B type, and one final R transformation. Let us give an explicit example. 
If (u~u + ) <G RG3, we will have 



Tsb{u 



RoBoG{u 



(9.12) 



where the explicit expressions of the R, B and G transformations have been given in Eqs. (|9.3|) . (|9.4j) . (|9.5j) above. 

Collecting all these ingredients together, it is easy to evaluate the probability Prx™ for a small-billiard era RX to 
contain a certain number n of epochs. These probabilities are obtained by integrating the w form over the suitable 
domain, as explained in Table IIXI Because the domains are not defined by straight lines, but by hyperbolas, the 
explicit expression for such probabilities would be somewhat awkward to write down. Nevertheless it is nevertheless 
always possible and straightforward to calculate it. 

As the most direct example, we evaluate the probability Prb 2 j i-e. the probability for a RB era to consist of 2 epochs. 
As we see from Fig. [S] and Table 11X1 the subregion of RB corresponding to this case is the simplified domain 



< u < (/), ul < u + < up, 

(j) < U~ < 1, U a < U + < Up, 



(9.13a) 
(9.13b) 



where u^ a = — (u +2)/(u +1), and u a = u^ a at u = <f>. The probability then reads 



P, 



RB 2 



du 



du+ 

(u+-u~y 2 



du 



(u+-u-) 2 



0.1240 



(9.14) 



We remark that, while, in the BKL case, the probability for an era to consist of k epochs was a monotonically 
decreasing function of fc, in the small billiard, this probability is a non-monotonic function of k: it increases for the 
lowest values of k and then starts to decrease. Furthermore, the probabilities Prb'G 71 for an era to contain n epochs 
in the RG U RB' domain can also be investigated. The most interesting information is obtained in the asymptotic 
region for large w + , as described in the above. For the RG region, in this asymptotic regime, the u~ boundaries can 
be considered as vertical straight lines, i.e. —1/2 < u~ < 0, and the u + boundaries can be considered as horizontal 
straight lines, as given by the asymptotic behavior of the functions u^ — > (n — l)/2, m" — > (n — 2)/2. On the other 
hand, the probabilities for the region RB' are obtained applying the transformation (|9.10[) to the appropriate functions 
in Table llXl For large values of u + , the region RB' is delimited by vertical straight lines u~ = — 1 and u~ = —1/2, 
and the horizontal lines delimiting the n-cpoch era starting region are u + = (n — 2)/2 and u + = (n — 3)/2. 
This way, the probability, say Prb'g™, to obtain a n-epoch era starting from the RG U RB' extended region is given 
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by the integral of the u form over the appropriate domain, i.e.: 



RB'G" 



±ln3 



n-l 
2 



n -2 
2 



diT 



du 



L/2 (U+-U ) 2 7"-^ 



n-2 
2 



-1/2 



(it+ - U~) 5 



For large values of n one finds 



4 1 

Prb'g™ — r~o^I- 
m 3 n z 



1 ^in-^- 1 ) 2 



In 3 n(n + 2) 

(9.15) 

(9.16) 



Note that the result (|9 . differs from the corresponding quotiented big-billiard (or BKL) result which is P m ~ j^^z- 

This difference shows that the natural notion of small-billiard era cannot be precisely mapped on the usual notion of 
BKL eras. 



B. Quotiented small billiard 



The direct, seemingly natural red- return map approach to the small billiard discussed in the previous subsection 
leads to a rather complex description of its dynamics. In particular, the fact that the n-epoch boxes RG n and 
RB n defined in the previous subsection have curved boundaries would make it difficult to define a small-billiard 
analog of the nice continued-fraction BKL description of the (big billiard) era dynamics. Here, we wish to show how, 
starting from the definition of the small-billiard, one can recover its deep relation with the big billiard, and uncover 
the technically simple BKL-like description of its dynamics. To start with, let us show several ways of relating the 
small-billiard dynamics to the big-billiard one. First, we can see this relation by using a graphical representation of 
the big-billiard dynamics which has been introduced by BKL. Namely, we mean the plot of the three logarithmic 
scale factors a = ft 1 = — In a, ft = ft 2 = — In 6, 7 = ft 3 = — lnc as functions of r. Modulo a conventional change 
of sign in the definition of a, /3, 7 (introduced here to ensure that a, /3, 7 are always > 0), this representation was, 
e.g., used in Fig. 2 of the review Q. In terms of this graphical representation, the dynamics of the variables 
a' , ft' , 7' of the small billiard is simply related to that of the variables a, ft, 7 of the big billiard in the following way: 
starting from a BKL graph of a(r), ft(r), 7(7-) [in which the three lines a(r), ft(r), 7(7-) can cross and keep (nearly) 
constant slopes except when one of them touches, from above, the horizontal axis (gravitation wall)] one can define the 
corresponding graph of c/(r), ft'(r), Y(t) simply by defining c/(r), for each r, as a' (r) = inf [a(r), ft(r), 7(7")], 7'(t) 
as 7'(t) = sup[a(r), /3(t), 7(7-)] and ft' (j) as the remaining middle curve among the three curves o>(t), ft(r), 7(7"). In 
other words, a'(r) is defined as the lower envelope of the three original a, ft, 7 curves, 7'(r) the upper envelope, and 
ft'(r) the intermediate curve. The three curves a'(r), ft'(r), 7'(r) satisfy a'(r) < ft'(r) < 7'(t) and keep constant 
slopes expect when either two of them "collide" (symmetry wall), or when the lower one, i.e. a'(r) "collides" with the 
a' = axis (a'-gravitational wall). One easily sees that a', ft', J 1 are the logarithmic scale factors of a small billiard 
(with gravitational wall a! = and symmetry walls a' = ft' and ft' = 7'). In other words, a suitable (time-dependent) 
re-ordering of the three big-billiard variables a, ft, 7 transforms the big-billiard dynamics in the small-billiard one. 
Reciprocally, starting from the graph giving the three curves o/(t), ft' (t), Y(t) of a small-billiard dynamics one can, 
starting from an arbitrary identification of a, ft,') with a',/3',7' at some initial time, extend the definition of the 
ct(r), ft{r), 7(7") curves by the condition that the only changes of slopes of these curves happen when one of them 
touches the zero axis. This graphical reasoning shows that there is an essential equivalence between the two billiards. 
Note in passing that this graphical approach can also clarify why the definition used in the previous subsection of the 
red-return eras introduces an artificial difference between the two billiards. Indeed, the red-return eras arc defined 
by the "collisions" ft'{r) = 7'(r) between the two upper curves. By contrast, the usual big-billiard eras are defined 
by a different condition consisting in finding when the "oscillations" between the crossing curves a(r) and ft(r) cease 
to give rise to oscillations between either a and 7 or ft and 7. The difference between these two definitions of "eras" 
give rise to the technical differences found in the previous subsection. 

A second way of relating the two dynamics consists of introducing some further quotienting of the small billiard. 
Explicitly, if we replace the small billiard by its kaleidoscopic version (as explained above in the big billiard case), 
i.e. by replacing the single moving ball of the small billiard by its six images with respect to the symmetry walls, we 
end up with a quotiented small billiard where the symmetry walls have no effect (because of the equivalence of the 
53 orbits) and where the only effective collisions take place on the gravitational walls. Finally, we conclude that the 
symmetry-wall-quoticnted small-billiard coincides with the symmetry-wall-quoticntcd big-billiard. This shows again 
that, modulo some discrete relabcllings, the two billiards are essentially identical. 
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X. BRIEF CONCLUDING REMARKS 



In this work, we have analyzed several aspects of "chaotic" cosmological billiards, in (3 + l)-dimensional gravity, and 
of their statistical behavior as one approaches the cosmological singularity. We have reviewed how the dynamics of 
the diagonal degrees of freedom (logarithmic scale factors /3 a ) of the spatial metric near the singularity is conveniently 
described by Lorentzian, or (after projection) hyperbolic-space billiards. We emphasized that the hyperbolic-space 
billiard table for the usual, Bianchi IX abc system is an ideal triangle (with three vanishing angles) , which contains six 
copies of the "small billiard" that is naturally "related" to the Weyl chamber of an hyperbolic Kac-Moody algebra, 
see Fig. [5] We reviewed several useful facts about integral invariants in Hamiltonian systems, and showed how their 
application to cosmological billiards allowed one to derive several forms and measures that are invariant under both the 
continuous and the discrete billiard dynamics. Contrary to previous treatments of cosmological billiards (starting with 
the classic work of Bclinski, Khalatnikov and Lifshitz, BKL), we did not use the six-fold symmetry group (S3) of the 
Bianchi IX a, 6, c system to symmetry-quotient its dynamics. This led us to defining a richer "hopscotch dynamics" 
between several sub- regions of the two-dimensional phase-space (u~,u + ) parametrizing successive Kasner epochs. 
Several aspects of this hopscotch dynamics have been discussed in detail: (i) the existence of a non-normalizable 
measure on the two-dimensional (u~,u + ) hopscotch court; (ii) the existence of a non-normalizable measure on the 
single variable u + , i.e. on the Kasner circle parametrizing the exponents of successive Kasner epochs; (iii) the existence 
of a normalizable measure on the subset of the hopscotch court describing the first epochs of successive eras; (iv) the 
link between the unquoticntcd hopscotch dynamics, and its quotiented version, equivalent to the usual BKL dynamics. 
Several statistical features of the hopscotch dynamics have been discussed, e.g. (1) the joint probability P ni ,n 2 f° r 
two successive eras to have specified lengths m and 712, and the fact that the random variables m and 712 are not 
statistically independent; (2) the "anisotropic" behavior of the hopscotch dynamics, i.e. the fact that the successive 
corners, between which the billiard ball representing the metric bounces, are statistically correlated. We briefly 
discussed the link between periodic orbits in the unquotiented hopscotch court, and the usually discussed periodic 
orbits in the quotiented, BKL description. Finally, we discussed the relation between the billiard dynamics within the 
full ideal triangle associated with the (diagonal) Bianchi IX model, and the dynamics between the "small billiard" 
that naturally arises in the treatment of the gravitational dynamics that uses an Iwasawa decomposition of the spatial 
metric. 



Appendix: Integral invariants in Hamiltonian systems 

As a complement to Section lTVl we recall here some (not always well known) facts about integral invariants in Hamil- 
tonian systems. When considering a general (possible time-dependent) Hamiltonian dynamics with (Hamiltonian) 
action (|4.1[) . one defines the Poincare-Cartan one- form 

a { pl := Pl dq l - H(q,p,t)dt, (A.l) 

defined in the (2n+ l)-dimensional extended phase space X^ 2n+1 ^ = {(q l ,pi, t)}. The Poincare-Cartan one-form (|A.1[) 
then defines a relative integral invariant of the unparametrized Hamiltonian flow. This means that the integral 

I{C)=j c a% (A.2) 

of en 1 ' over any closed curve C in extended phase space remains the same if one displaces C in an arbitrary manner 
along the unparametrized flow lines of the Hamiltonian dynamics. [By unparametrized flow lines we mean here the 
unparametrized lines in X^ 2n+lS> whose tangents are everywhere parallel to the Hamiltonian flow, i.e. proportional 
(without being necessarily equal) to {dH/dpi, —dH/dq 1 , 1).] An equivalent formulation of this property is to say that 
the Poincare-Cartan two-form 

: = -do-^l = dq i A d Pl - dt A dH(q, p, t) ( A.3) 

is an absolute integral invariant of the unparametrized Hamiltonian flow, i.e. that the integral of oj( 2 > over an arbitrary 
two-surface £ (with boundary) is invariant as £ is moved, in an arbitrary manner, along the unparametrized flow 
lines. Note that these invariance properties are stronger than those that are usually stated, which only refer to the 
invariance of the simpler 'Liouville' forms in the unextended phase space X^ 2n ^ = {(q l ,pi)}, 

4 1 ' ~ V,dq\ (A.4a) 
w[ 2) := -da^ = dq l A dpi, (A.4b) 
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under the time-parametrized Hamiltonian flow (dH/dp, —dH/dq, 1). 

Moreover, the absolute invariance of the two-form u>pq, Eq. (|A.3[) . allows one to construct an invariant measure on 
any (2n)-dimcnsional transverse section moving along the unparametrizcd Hamiltonian flow, namely 



Opp := c(n) ^ppj = c(n)uipc A ... A Wpp (with n factors). (A. 5) 

Here c(n) is a numerical factor which can be taken to be (-^("-^(n!)" 1 if one wishes to recover the usual normal- 
ization. In the case when one restricts oneself to displacements along the time-parametrized Hamiltonian flow one 
can replace the extended phase space form Slp^Q by the Liouville measure 

VL { l n) := c(n) (w[ 2) ) A " = <V A d 1 2 A - A d 1 n A d Pi A d P2 A - A d Pn (A.6) 
on the unextended phase space X^ 2 ^ = {(g 1 ,^)}. 



In addition, if one is considering a time-independent Hamiltonian H(p, q), and if one wishes to restrict oneself to the 
dynamics on a specific (2n — l)-dimensional energy hypersurface, say £ E , satisfying H(p,q) = E in unextended 
phase space X^ 2n \ the above results simplify in that one can drop the i7-dependent contribution in (|A.1[) and (|A.3[) 
(because dH(q,p) = on the energy shell) and conclude that the simpler Liouville-type two-form Eq. (|A.4bp . 
is invariant not only under the usual time-parametrized Hamiltonian flows q = d p H, p = —d q H in £^ 2n ^ 1 \ but also 
under more general 'many fingered time flows' q = Fd p H , p = —Fd q H, involving an arbitrarily varying time-rescaling 
function F(q,p,t). 

This leads to introducing two possible measures associated to the dynamics on the energy-surface £ E n 1 \ A first 
measure is the standard energy-shell reduced Liouville measure 

^ 1] =^l n) 6m<l,p)-E), (A.7) 

which is a {2n — l)-form, and yields a smooth measure on the (2n — I)-dimensional energy shell £( 2n_1 ). 



A second possible construction (which is linked to the general theory of the reduction of phase spaces with symmetry) 

L , Eq. (|A.4b[) . under arbitrary 'glidings' along 



is to use the invariance of the two-form Wpp, Eq. (|A.3[) . or simply a/ 2 ^ 



the Hamiltonian flow (which takes place within £^ 2n ^) to define both a reduced (symplcctic) 2-form 



J red 



,( 2 ) 



(A. 



and a corresponding measure 



O 



(2n-2) 
red 



= c(n - f ) 



,(2) 



A(n-l) 



(A.9) 



(2n-2) _ 



£ E jTn (where Fh denotes the unparametrized Hamiltonian 



on the (2rt — 2)-dimensional quotient space Q E 

flow on fi 2n 2 ^). In other words, £ E n /J^h is the space of unparametrized Hamiltonian motions on £ E n x \ 



A 



concrete representation of this quotient space can be obtained by considering any transverse section of Fh on £ ( 2n 1 ', 
i.e. any 'initial conditions' for Fh- Note that this transverse section does not need to be taken at some fixed time t, 

■ (on £( 2n_1 )) under Fh then 



(2) 

but can have an arbitrary 'slope' in extended phase space. The invariance of ui P q 



guarantees that the (2n — 2)-form (|A.8p lifts to the quotient space Q 
the (2n — I)-dimensional space of (unparametrized) motions with energy E. 



(2n-2) _ £,(271-1) 



E 



E 



/Fh, i-e. defines a measure on 



Let us mention that there is a simple link between the reduced (2n — 2)-form SI 



(2n-2) 



Eq. (|A.8[) and the usually 



considered energy-shell reduced Liouville measure fl E n E , Eq. (|A.7|) . First, note that these forms define measures on 



different spaces: f^ 2 ^ 1 - ) 'lives' on the (2n — l)-dimensional energy surface £ E 



(2n-l) 



while Q 



(2n-2) 



1 (2n-l) 



of £§ n by the Hamiltonian flow Fh- To see the link between these constructs, we 



dimensional quotient Q E ui <^ E 
can introduce (at least locally) in the full, ambient (2n)-dimensional phase space X^ 2n 



'lives' on the (2n — 2)- 

these constructs, we 
{(g l ,Pi)} a new canonical 
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coordinate system where the n-th momentum coordinate p" cw is equal to the Hamiltonian: H(p° ld , q id) — p" cw - in 
this new canonical coordinate system, the n-th conjugate position coordinate <7™ cw is such that 

= dH = 1 (A 10) 

dt dp™ w ' v ■ J 

while the remaining coordinates (q z ,pj), with i = l,...,n— 1 all satisfy q l ncw = = pf cw , i.e. are invariant under the 
Hamiltonian flow Tn (considered in unparametrized form, i.e. with many fingered time displacements: At = F(p, q)). 



If we denote the conjugate pair (9," ow ,p" cw ) simply by (s,H), and the (n — l)-other pairs by (<f,Pi), we see that 
the symplectic form ui^, Eq. (|A.4b[) . in the ambient phase space X< 2 ") reads 

4 2) (q, p) = oj£l(q, p) + ds A dH, (A.ll) 

where 

n-l 

1=1 

is clearly equal to the reduced symplectic form (IA.8I) on the quotient space Q^"~ 2 ' = £% n ~ jTu (independently of 
the values of E and s). 



If we now insert (|A.11[) in the general definition (|A.7[) of the energy-shell-reduced Liouville measure, one easily sees, 
using S(H - E)dH = 1, that 

^T ] =^tT 2) ^d S , (A.13) 

where, as explained above, s is a phase-space coordinate which is canonically conjugate to the Hamiltonian (which 
implies that ds/dt = 1 along the Hamiltonian flow) . 



So far, we have been considering any autonomous Hamiltonian system (dH/dt = 0). We can, in particular, apply 
the above results to general billiard systems, i.e. to a Hamiltonian of the form 

n 1 

H(p,q) = + VM, (A.14) 

where g tJ (q) is the matrix inverse of some (pseudo-)Ricmannian metric gij(q)dq'dq J , and where the (formal) potential 
function V oa (q) is equal to zero within some domain, say B (the "billiard table"), of the q variables, and equal to +oo 

(2) 

outside of this domain. In that case, the general invariancc of the reduced symplectic form under arbitrary, many 
fingered time motions on the energy hypersurface (see Eq. (|A.lip ) can be concretely used to show that the restriction 

(2) 

of the ambient phase space symplectic form uj l (q,p) on the boundary dB (with the constraint H(q,p) — E) of the 

billiard , say w rC s tr ('? restr :?'rcstr), is invariant both under each collision on dB and under each 'free flight' between 
two successive collisions. In other words, we are here considering transverse (Poincare-type) cross sections of the 

energy-reduced Hamiltonian flow on £^ n 1 \ defined, in a 'stroboscopic' manner, by the successive collisions. This 
allows one to extract from the continuous Hamiltonian flow (f> t [x(t) = 4>t(x(0)) with x = (q,p)j, the discrete 'billiard 
map' say T, such that xn+i = T(xn) where xn = (qN,PN) is the phase-space position just after the iV-th collision 
on dB and T(x) = 4> T ( x )+o(x) the stroboscopic Hamiltonian evolutional between two successive collisions (including 

the 'reflection' effect of the second collision). In other words, the restrictions oj^l t1 .(q Tes tr, p-rostt) of w^ 2 ' on the (2n — 2)- 
dimensional phase space of [dB] H ^ p q )—E after each collision gives us an infinite collection of concrete realizations of 

the reduced 2-form on the abstract quotient space Q^ 2n ^ = £^ 2n+1 ^ /Tr ■ It also yields several absolute integral 

(2) 

invariants of the billiard map T, namely w 1 . 0Str ('?restr,Prestr) itself and all its exterior powers, and notably the reduced 



Note that the time t(x) between two successive collisions generally depends on the starting position x, so that we need to use here the 
invariancc of the (Poincarc-Cartan J A.3t or Liouville llA.4bll ') two- form under the unparametrized many-fingcrcd Hamiltonian flow. 
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measure f2 
measure fl 
e.g., M )• 



(2ra-2) 
red ' 
(2n-2) 
red 



Eq. (|A.9[) . Note that the link (|A.13I) between the energy-shell Liouvillc measure and the reduced 
(invariant under the billiard map T) is well-known in the mathematical literature on billiards (see, 



Let us now discuss the various ways in which the above results can be applied to cosmological billiards. We start 
by recalling that the dynamics of the 'diagonal' degrees of freedom (i.e. the logarithmic scale factors f3 a entering the 
Iwasawa decomposition (|2.4p ) is described, near a cosmological singularity, by an action of the general form plj 



S B = / dx° 



where $ a = 4^0, and where V(/3) is a sum of 'exponential walls' 



NV{p) 



]Tca cxp {-2w A {(3)) 1 



(A.15) 



(A.16) 



with linear forms w A (/3), see Eqs. ([2.121 |2~T5I I2.16[) . The spatial gradients of the metric and of the other fields enter 
only in the coefficient ca of the exponential walls. In the near-spacelike-singularity limit ('BKL limit') the time- and 
space- dependent coefficients ca{x°, x) tend to some finite limits so that one can describe the asymptotic dynamics of 
the f3 a (x°, x) at each point of space by means of the Toda-like billiard (|A.15[[A".16[) (with ca replaced by their limits). 
A further approximation (which also holds in the BKL limit) consists in replacing the exponential walls (|A.16[) by 
their formal 'sharp wall limit', namely 



V O0 (/3)=J2^oo(-2waW)) 



(A.17) 



where the formal sharp- wall Ooo-function is defined as: Ooo(x) := if x < and &oo(x) '■= +oo if x > 0. 

The action (|A.15|) with V(f5) — > Voo(f3) given by (|A.17[) defines a Lorentzian billiard dynamics in the /3-space. This 

dynamics can be equivalently described by the Hamiltonian action 



Sn 



dx° nJ a -H (p a ,7r a ) , 



(A.18) 



iG ab 7r Q 7r b + y oo (/3) 



(A.19) 



which is of the general type (|A.14|) (with a flat Lorentzian-signature metric G a b)- Here, n a denotes the conjugate 
momentum of j3 a , i.e. 



1 

N ab dx°- 



(A.20) 



Note that ir a is invariant under the redefinitions of the time coordinate a; (which affect both dx° and TV = N/^/g 
but leave invariant the product Ndx° tx Ndx°). 



As the (rescaled) lapse N is a Lagrange multiplier in the action Sp, Eq. (|A.15|) . we have the well-known Hamiltonian 
constraint stating that Hp must vanish, i.e. that we must constrain ourselves to the specific energy hypersurface 



Hp(p,Tr)=Ef, = 0. 



(A.21) 



We are therefore in the condition where we can apply the general results recalled above. [For simplicity, we can 
assume that we are working in any gauge where N is given as some autonomous function of /3 and tt. This is the case 
both of the r-time gauge N = 1, Eq. (O, and of the T-timc one N = p 2 = -G ab (i a (i b , Eq. (j2~25j) . ] In particular, 
we see that 

uf ] = -d[TT a d(3 a ] = df3 a A dn a (A.22) 
is an absolute integral invariant of the Hamiltonian flow, as well as the corresponding energy-shell measure 



(2d- 
E«=0 



, 1} =c(d) (u:fy d 5{H p ) 



(A.23) 
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In addition, we can consider the (double) reduction of the 2-form aA on the quotient space Q\ 



o Ef,=0 /•'Hp, 



ufjed = W 1 A dTr a } Q (2d-2) , 



(A.24) 



and the corresponding measure on Q^ d 2 ' 



(2d-2) 
/3red 



c(d-l) 



/a) 



A(d-l) 



(A.25) 



The latter reduced measure is related to the Liouville-type measure (|A.23[) via the general result (|A.13[) where s is a 
/3-phase-space coordinate which is canonically conjugate to Hp, Eq. ()A.19j) . These results are particularly simple if 
one uses the r-time gauge where N = 1, so that the Hamiltonian (|A.19| is the sum of a constant-Lorentzian- metric 
kinetic term iG af> 7r a 7r;, and of a sharp-wall billiard potential. 

The integral invariants we have just discussed concern the dynamics of the Lorentzian billiard in /3-space. They 
would be useful if one were studying the full /3-space billiard dynamics. However, in this paper, we are interested 
in discussing the projection of the /3-billiard on the hyperbolic space W4-1, e.g. described by the unit hypcrboloid 
(|2.26|) in /3-space. This projection is obtained by separating out the motion along the 'radial direction' p, Eq. (|2.23[) . 
Indeed, setting N = N / p 2 and A = lap, the billiard action (|A.15[) , (|A.17|) can be rewritten as (see, e.g. [l4| for the 
d = 3 case and [13] for the general case) 



S 



1 

2N 



-X 2 + G ab n b )-NV 00 (j) 



(A.26) 



the 



where Voo(j) = J2a ®oo(— ^wa(j))- [As usual, we are using here the simplifying fact that, in the limit p — > - 
potential term becomes independent of p, i.e. becomes T-time independent.] 

In the T gauge (|2.25|) . i.e. N = N / p 2 = 1, the radial kinetic energy term — |A 2 decouples from the angular motion 
terms and leads to a uniform radial motion: dX/dT = const, in that gauge, one can simply work with the "angular 
action" 



dT 



1 d^d^ 
2 Gab lT ~dT _Ko(7) 



submitted to the constraint that the (constant) angular-motion energy 



1 d~f a d"f b 



(A.27) 



(A.28) 



be equal to (dX/dT) 2 = const. 



Note that while the /3-space action Sp, Eq. (|A.18[) , corresponded to a phase-space (/3 a ,7r a ) with 2d dimensions, 
the reduced action S 7 , Eq. (|A.27|) . corresponds to a phase space with 2{d — 1) dimensions. In order to explicitly 
describe the reduced dynamics, one needs to choose some parametrization of the {d — 1) dimensional hyperbolic space, 
say q l , where the index i takes only d — 1 values. The hyperbolic metric on Hd-i will have some expression, say 
= l,...,d- 1) 



G ab d7 a dj b 



ds 2 



g l:i {q)dq l dq J , 



and the angular action (|A.25[) will read 



dT 



1 dq l dq> 

2 9ij[q) dT^T 



Voo(l(q)) 



The conjugate momenta to the q l, s read 

while the angular-motion Hamiltonian will read 

H~ t (q\p l ) 



Pi = 9ij(l) 



dq j 
~dT 



-,9 ij (q)PiPj + Voo(l(q)), 



(A.29) 



(A.30) 



(A.31) 



(A.32) 
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where g l -'{q) denotes the inverse of the (covariant) metric gij(q)- 

Similarly to the discussion above, we can now introduce the (reduced) Poincare-Cartan one-form 



and the corresponding two-form 



a' 1 ' := Pl dq l - H 7 (q, p)dT, (A.33) 



= -da {1) = dq l A dpi -dTA dH^(q,p). (A.34) 



A 2) 



As before (respectively u)^) defines a relative (resp. absolute) integral invariant of the unparametrized Hamilto- 

(2) 

nian flow in extended phase-space (q,p,T). And, as before, we can use the absolute invariance of the two- form w 7 
to construct an invariant measure. As we are again in a situation where we can work on a fixed-energy hypcrsurfacc 
(here H 7 = E~, = const, after eliminating the uniform radial motion A(T)) we can drop the last, independent term 

in io^\ Eq. (|A.34|) . and work with the usual (7-spacc) symplcctic form w 7 2 ' = dq l A dpi. As before we end up with 

having a whole set of integral invariants of the billiard dynamics in 7-space (i.e. on the hyperboloid Hd-i)' the 
(2) 

two-form uij itself, and its various exterior powers, and its energy-shell restricted measure 



(2d-3) _ 



c(d - 1) (^ 2) ) A<d ^ d(H 7 (p,q)-E 7 ). (A.35) 

3r, we can also use the invariance o 



Moreover, we can also use the invariance of the reduction of the symplcctic form on the quotient space Q^ d ^ 



(2) 

UJ, 



7red 



[dq l A dpi] Q (2d-4) (A.36) 



and its maximal exterior power 



^=c(d-2)(JZ) . (A.37) 



In addition, we still have a link of the type (|A.13|) (where s is a 7-phase-space coordinate such that ds/dT = 1 
along the Hamiltonian flow), and we also know that the abstract quotient-space reduced symplectic form oj^ red can 
be concretely computed by restricting dq l A dpi by two conditions: H y (q,p) = _E 7 and any cross-section condition 
transverse to the Hamiltonian flow. In particular, we can use the events of collisions on successive walls of the billiard 
as cross-sections, and thereby prove that w^ onision and ^collision are invariants of the discrete 7-billiard map T which 
connects a collision to the next. 



In Section (|IVp. we apply these general results to the case of the BKL cosmological billiards, in d = 3 spatial 
dimensions, working within the radially-projected picture on the 7-space, i.e. on the hyperbolic plane 7^2- 
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